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Abstract 



We show how the Hintermann-Merlini-Baxter-Wu model (which is a gen- 
eralization of the well-known Baxter- Wu model to a general Eulerian triangu- 
lation) can be mapped onto a particular infinite-coupling-limit of the Ashkin- 
Tcller model. We work out some mappings among these models, also including 
the standard and mixed Ashkin-Teller models. Finally, we compute the phase 
diagram of the infinite-coupling-limit Ashkin-Teller model on the square, tri- 
angular, hexagonal, and kagome lattices. 

Key Words: Baxter- Wu model; Ashkin-Teller model; Hintermann-Merlini-Baxter- 
Wu model; Infinite-coupling-limit Ashkin-Teller model; partial trace transformation; 
plane Eulerian triangulation. 
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1 Introduction 



In 1944 Onsager P] solved the nearest-neighbor square-lattice Ising model without 
magnetic field. This constitutes a major milestone in Statistical Mechanics. This 
achievement opened the door to find other exactly two-dimensional solvable models 
using different approaches. Two of such lines of research were (1) to consider models 
with multi-spin interactions rather than nearest-neighbor couplings, and (2) to look 
for models with larger symmetry groups than Z 2 . 

Two different groups obtained positive results along the first line of research using 
3-spin interactions. In particular, in 1972 Hintermann and Merlini [2] solved the Ising 
model with 3-spin interactions on the union-jack lattice. They use a mapping of this 
model onto the 8- vertex model [3J, and used Baxter's solution of the latter [3] (another 
milestone in Statistical Mechanics) to derive the critical points of the model. Soon 
after this result, Baxter and Wu |U[S] found in 1973 the exact solution of an Ising 
model on a triangular lattice also with pure 3-spin interactions. The method chosen 
was the Bethe Ansatz [3]. This model is very interesting, as it is believed [6] to belong 
to the same universality class as the 4-state Potts model [3j[7], but it does not show 
the logarithmic corrections displayed by the latter model [8TfTU] . The field theoretical 
interpretation of these differing behaviors is the following. The dominant part of 
the action describing the continuum limit of both the 4-state Potts model and the 
Baxter- Wu model is that of a free bosonic field. Both models contain further strictly 
irrelevant operators, which are not necessarily identical, since they may depend on the 
discrete symmetries of the spins and of the lattices on which the models are defined. 
But crucially, the 4-state Potts model contains in addition a marginally irrelevant 
operator, which is not present in the Baxter- Wu model. This latter operator induces 
logarithmic corrections under the renormalization-group flow that takes the 4-state 
Potts model to the fixed-point theory. Alternatively, one may say that the Baxter- Wu 
model is designed so that the amplitude of the marginally irrelevant operator is set to 
zero. The Baxter- Wu universality class is also described via its universal amplitude 
ratios in Ref. |llj . 

Both models, the Hintermann-Merlini (HM) and the Baxter-Wu (BW) models, 
had many features in common: they were defined on plane triangulations (i.e., all 
faces were triangles), the dynamical variables were Ising spins, and the Hamiltonian 
only included 3-spin terms (one for each triangular face). These two models belong 
to a more general family of models, that we will denote as the Hintermann-Merlini- 
Baxter-Wu (HMBW) model. We define this model on any plane triangulation G = 
(V, E) of vertex set V and edge set E. Then on each vertex x G V we place an Ising 
spin a x = ±1. These spins interact through the Hamiltonianfj] 

WhmBW = —J <Ji&j&ki (1-1) 

A={i,j,k} 

where the sum is over all triangular faces A = {i,j, k} of G, and J is a coupling 

1 We could define the HMBW model with face-dependent couplings; but we will not need this 
generalization in this paper. This comment also applies to the other models considered here. 
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constant. This kind of models (and generalizations of them) were already considered 
by other authors [T2|[T3] from a group-theoretical point of view. In this paper, we 
will focus on the particular case of Eulerian triangulations: a graph is Eulerian if the 
degree of all its vertices is even. This property will play an important role in the 
following sections. 

The second research line produced (among others) the so-called Ashkin-Teller 
(AT) model [14J . This model was introduced in 1943 and generalized the Ising model 
to a 4-state model. It was soon recognized that it could be written as two coupled 
Ising models [T5]. More precisely, given a graph G = (V, E), we place on each vertex 
x G V of the graph two spins cr*, 7$. These spins interact via the Hamiltonian: 



where the sum is over all edges of the graph. We are not aware of any exact solution 
to this model, or to any of its particular symmetric case, to be reviewed in Section [3]). 
This AT model contains, as particular cases, the Ising and the 4-state Potts models. 

Furthermore, in the symmetric case K' 2 = K2, there is a curve in the (K2, K4)- 
plane that can be mapped to a soluble 6- vertex model [3]. This curve is self-dual, and 
part of it is critical, with critical exponents varying continuously as we move along 
that curve. This phenomenon provides the simplest counter-example to the usual 
notion of universality. 

At first sight, the relation between the HMBW and the AT models is rather weak: 
both models can be defined in terms of Ising spins, and we can find the 8-vertex 
model in both the solution of the HM model, and in the computation of the self-dual 
curve of the symmetric AT model. One of the goals of this paper is to show that 
the relation is deeper: there is an exact mapping between the HMBW model and an 
infinite-coupling-limit of the AT model (ICLAT), to be defined below. In particular, 
certain curves in the phase diagram of the ICLAT model on a graph G are equivalent 
via those mappings to a HMBW model on a certain triangulation G'. 

In order to achieve our main goal, we need to work out a series of exact mappings 
among the HMBW, ICLAT, AT models and some other models we will introduce 
in Section |3J Some of them are already known results in the literature, others are 
alternative versions of known results, and finally, some of them are new. 

As an application, our last goal is to compute the phase diagram of the ICLAT 
model on the square, triangular, hexagonal, and kagome lattices, which include one or 
two HMBW models as particular cases. In all cases, we find a point which corresponds 
to the point where the self-dual curve of the AT model on the corresponding lattice 
hits the ICLAT plane. This point is denoted B in Section |5j For the first three 
lattices, point B belongs to a subspace of the ICLAT model that can be mapped 
onto the HMBW model. So we can approach it following two independent paths: 
either along the HMBW subspace, or by considering the limit of large couplings of 
the self-dual curve of the corresponding AT model. For these lattices we can make 
predictions of the critical exponents along these two paths, and it turns out that they 
differ. 




(1.2) 



e={ij)£E 
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The paper is organized as follows. In Section [2] we will show the graph-theoretic 
set-up. In Section [3] we will define carefully the distinct models we will consider in 
this paper. In Section 2] we will work out the mappings between the AT and mixed 
AT models, the mapping between the mixed AT and the ICLAT models, and finally 
the partial trace transformation. This latter mapping allows us to relate the HMBW 
and the ICLAT models on certain lattices. We conclude this section with a nice 
mapping between the BW model and a loop model. The next Section [5] is devoted 
to the description of the ICLAT phase diagrams for several common lattices: square, 
triangular, hexagonal, and kagome. In addition we provide the critical exponents 
for the limiting point of the AT self-dual curve when the couplings tend to infinity. 
Finally, in Section [6] we summarize our findings. In Appendix |A] we describe the 
decimation transformation and how it is related to the partial-trace transformation. 

2 Graph-theoretic preliminaries 

In this section we will introduce the main lattices where we will define our physical 
models of Section [3j 

Let us consider a graph G = (V, E) with vertex set V, edge set E, and embedded 
in the plane (i.e, G is planar). We can always define the dual graph G* = (V*,E*) 
of G in the standard way as follows: To each face / in G, there corresponds a dual 
vertex /* G V*; and for every edge e G E, we draw a dual edge e* G E*. If the 
original edge e lies on the intersection of two faces / and h (possibly / = h), then 
the corresponding dual edge e* joins the dual vertices f*,h* G V*. We can draw G 
and G* in the plane in such a way that each edge e G E intersects its corresponding 
dual edge e* G E* exactly once. 

A graph G is a plane triangulation if it is planar and if every face of G (including 
the outer one) is bounded by a triangle. If G is a plane triangulation, then G* is a 
planar cubic graph (i.e., a planar 3-regular graph). Therefore, if G is an Eulerian 
plane triangulation (i.e., all vertices have even degree), then G* is a bipartite cubic 
graph. This implies the following theorem due to Heawood [161 Exercise 18, Chapter 6] 
(see also [17]): 

Theorem 2.1 A plane triangulation is 3-colorable if and only if it is Eulerian. 

Remarks. 1. In the above theorem, 3-colorable means that we can color the vertices 
of G with three colors in such a way that two adjacent vertices are not colored alike. 

2. This theorem holds for arbitrary Eulerian plane triangulations. In particular, 
we allow for multiple edges. Let us show a simple example. For any k > 2, the graph 
Gk consists of a pair of vertices connected by 2k paths, alternating (as one draws 
the graph on a plane) of lengths 1 and 2. The smallest example Gi is depicted in 
Figure [TJ 

3. Given a Eulerian plane triangulation, then it is uniquely 3-colorable (modulo 
global color permutations). 
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Figure 1: Graph G2 which is an Eulerian plane triangulation. 



As a result, given any Eulerian plane triangulation, there is a unique way to split 
the vertex set into three subsets V = V\ U V2 U V3, such that they are mutually 
disjoint (i.e., Vi fl Vj = for all j 7^ i)), and for each edge e G E, its endpoints 
belong to different subsets (i.e., if e = {a, b} G E, then a G Vi and b G Vj for 
i ^ j). Furthermore, the edge set can also be partitioned into three disjoints sets 
E = E\2 U -E13 U E23, such that any edge e = {a, b} G E^ has one of its endpoints in 
Vi and the other one in Vj with i ^ j. 

The dual G* = (V*,E*) of an Eulerian plane triangulation G = (V,E) is a 
bipartite cubic graph. Therefore, we can properly color the vertices of the dual graph 
G* with two colors, or equivalently, we can color the faces of G with two colors in 
such a way that any two neighboring faces (i.e., two faces with a common edge) are 
not colored alike. Therefore, this proper face 2-coloring induces a partition of the 
(triangular) faces F = F(G) into two disjoint sets F = F\ U F 2 , such that all all 
faces in Fi are colored i. For a plane triangulation G, the following identity follows 
from the handshake lemma on G*: \E\ = 3|F|/2, which implies using Euler's formula 
l-^l - 1^1 = 1-^1 — 2, that the number of faces \F\ is even: \F\ =2(\V\— 2). Furthermore, 
\Ei\ = 1-^2 1 = \V\ — 2, because G* is 3-regular. 

A graph G is a plane quasi-triangulation if it is planar and if every face of G 
(except the outer one) is bounded by a triangle. This outer face is bounded by a 
cycle of length i > 4. Plane quasi-triangulation are useful when considering a finite 
piece of a regular lattice with free boundary conditions: see e.g., Figure [2J If G is a 
plane quasi-triangulation with an outer face bounded by a cycle of length I > 4, we 
have to make some modifications to the above results. On one side, if the degree of 
all inner vertices is even, then the quasi-triangulation is 3-colorable. On the other 
hand, if there are no bridges it is also uniquely 3-colorable (modulo global color 
permutations). 

The number of faces is \F\ = |.F in | + 1, where |F in | is the number of inner triangular 
faces. Then the number of edges is given by l^l = (3|F in | +£)\/2. Therefore, |F in | and 
£ should have the same parity, as \E\ is an integer. The number of triangular faces 
is given by |F; n | = 2\V\ — i — 2. Therefore, we need an outer cycle of even length to 
have an even number of inner triangular faces. 

The dual of a plane quasi-triangulation with all the inner vertices of even degree 
is a graph G* = (V*,E*) constructed as follows: the vertex set V* = V^ U {^outl 
contains a vertex v* G V|* for each inner triangular face in G, and an additional 
vertex v* nt accounting for the outer face in G. The former vertices have degree three, 
and the latter has degree i. The edge set E* is constructed in the usual way. The 
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subgraph induced by the inner vertices V* n is indeed bipartite, so the set of inner faces 
F in can also be split into two disjoint sets F in = F\ U F 2 . However, we cannot ensure 
that both sets Fi have the same cardinality, as the dual graph G* is not regular. 

We are going to study two classes of plane triangulations: the standard trian- 
gular lattice (which is a 6-regular graph), and the infinite family of Eulerian plane 
triangulations introduced in [T8] . 

2.1 Triangular lattice 




Figure 2: Triangular lattice of size 7x7 with free boundary conditions. The vertex set 
has a partition into three disjoint sets, depicted as gray, white, and black dots. The 
edge set is formed by the black and gray edges. The black edges show the subgraph 
G\i = (Vi U V 2 , E 12 ), which is a hexagonal lattice. 

An infinite triangular lattice G = (V, E) is a 6-regular triangulation of the infinite 
plane, as shown in Figure [21 It is uniquely 3-colorable, therefore we can split the 
vertex set into three disjoint subsets: V — V\ U V 2 U V3. Note that for each % 7^ j, the 
subgraph = (V, U Vj, Eij) forms a hexagonal lattice. 

The dual lattice G|- = (Vk,E*) contains as a vertex set the third subset (k 7^ 
i, j) of V. This also works on the other way around: for each i — 1,2, 3, the vertex set 
Vi can be regarded as the vertex set of another triangular lattice (with larger lattice 
spacing), and whose dual is precisely the hexagonal lattice Gjk = (Vj U Ejk) with 
i 7^ j, k. For instance, in Figure [2] the white and gray dots and all the black edges 
form a hexagonal lattice, whose dual has a vertex set formed by the black dots. 

A finite piece of the triangular lattice is a finite subset of the infinite lattice 
described above with linear dimensions L x x L y and free boundary conditions. Figured 
shows an example with L x = L y = 7. This is a plane quasi-triangulation; but it is 
not Eulerian, as there are precisely two vertices with odd degree (the top right and 
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the bottom left corners). All the inner vertices % have degree d{ = 6, but the degree 
of the vertices belonging to the outer cycle can be 2, 3, or 4. The length of this outer 
cycle is £ = 2(L X + L y — 2). In order to obtain an Eulerian plane triangulation from 
this graph, we proceed as follows: 

(a) Add an edge / joining the two vertices with odd degree. We obtain an Eulerian 
graph with all faces being triangles, except two, which are bounded by cycles 
of length £' — L x + L y — 1. 

(b) If £' is not a multiple of 3, then subdivide the extra edge with 1 (resp. 2) 
vertices when £' = 2 (mod 3) (resp. £' = 1 (mod 3)). So we end with two faces 
bounded by cycles of length multiple of 3, and all vertices along this cycle have 
even degree. 

(c) Take one cycle, and number the vertices 1,2,3 according to the sub lattice Vi 
the vertex belongs to. Then form triangular faces by adding one edge between 
two vertices labeled 1 and 3, with a single vertex labeled 2 in between. 

(d) We obtain one cycle of even length with all vertices of odd degree (as in the 
previous step we added a single edge to each vertex). The vertices on this cycle 
are labeled alternatively 1 and 3. Now place a new vertex (labeled 2) inside the 
face bounded by this cycle. Then add edges between this new vertex and any 
of the vertices defining the cycle. Indeed, the degree of all vertices is now even. 

(e) Repeat the last two steps on the other cycle of length multiple of 3. 

In this way, from a finite subset of linear size L x x L y of a triangular lattice with 
free boundary conditions, we have produced an Eulerian plane triangulation. Indeed, 
the added vertices and edges are a negligible fraction of the total number of vertices 
and edges, when L x ,L y — > oo. Therefore, we expect that their contribution to the 
infinite-volume free energy would be zero. 

2.2 A family of Eulerian plane triangulations 

In this section, we follow Ref . [18] . Let us construct this family of Eulerian planar 
triangulation by starting from a connected planar graph G = (V,E). We can then 
define its dual G* = (V*, E*). The next step is to build the graph G = (V U V*, E) 
with vertex set KUF* and edges ij whenever i G V lies on the boundary of the face 
of G that contains j G V*. The graph G is a plane quadr angulation: on each face 
of G, one pair of diametrically opposite vertices corresponds to an edge e G E, and 
the other pair corresponds to an edge e* G E*. In other words, G is the dual of the 
medial graph M.(G) = M.(G*). Conversely, every quadrangulation G arises via this 
construction from some pair G, G*. (See [T8J Figure 1].) The vertices of G (resp. G*) 
are depicted gray (resp. white), and the edges in E are depicted as solid black lines 
in Figure [31 

Let us now define the graph G constructed from G by adjoining a new vertex in 
each face of G, and four new edges connecting this new vertex to the four corners of the 
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Figure 3: Two pieces of planar triangulations with free boundary conditions: (a) 
the union-jack lattice of size 5x5, and (b) the bisected hexagonal lattice of size 
4x4 (measured in units of the underlying triangular Bravais lattice i.e., the gray 
dots). The vertex set has a partition into three disjoint sets, depicted as gray, white, 
and black dots. The gray and white dots and the black solid lines form the lattice 
G: square (a) and diced (b). The new vertices and edges added to G to get G are 
depicted as black dots and solid gray lines, respectively. 



face. The new vertex is depicted as a black dot, and the new edges as solid gray lines 
in Figure EJ The graph G is an Eulerian plane triangulation, with vertex tripartition 
V = VUV*UV 3 , where V 3 is the set consisting of the new degree-4 vertices. Conversely, 
every Eulerian plane triangulation in which one sublattice consists of degree-4 vertices 
arises in this way. Indeed, this does not cover the triangular lattice. 

If G = G* = G = square lattice, then G = union-jack lattice, displayed in 
Figure Eta). If G = triangular lattice, G* = hexagonal lattice, then G = diced 
lattice, and G = bisected hexagonal lattice, displayed in Figure Mjo)- According to 
the standard notation [19], the union-jack lattice is the [4 • 8 2 ] Leaves tiling, and the 
bisected hexagonal lattice corresponds to the Leaves tiling [4-6-12]. 

The infinite union-jack lattice is not regular since di = 8 for all i G V U V*, while 
di = 4 for all i G V3. This lattice can be regarded as a square Bravais lattice with a 
two-point basis. Notice that the union-jack lattice satisfies the following properties: 

• The graph G = (V U V*, E) has a square-lattice dual G* = (V 3 , E*). 

• The subgraph G 13 = (VUV3, £13), where E i3 are the set of edges joining vertices 
e = {ij} with i G V and j G V3, can be regarded as a decorated square lattice 
tilted 45° with respect the original lattice. Its dual is another tilted square 
lattice G{ 3 = (V*,E^ 3 ), where every edge in G{ 3 is doubled. The same occurs 
to the subgraph G 23 = (V* U V 3 , E 23 ). 



9 



The bisected hexagonal lattice is also non-regular: the degrees of the vertices in 
V (gray dots), V* (white dots), and V 3 (black dots) are respectively 12, 6, and 4 (see 
Figure 13(b)). This lattice can be regarded as a triangular Bravais lattice (formed by 
the vertices in V) with a 6-point basis. Notice that the bisected hexagonal lattice 
satisfies the following properties: 

• The graph G = (V U V*,E) is a diced lattice, and its dual is G* = (V 3 ,E*) 
which is a kagome lattice. 

• The subgraph G13 = (yuV3, -E13), where Ei 3 are the set of edges joining vertices 
e = {ij} with i G V and j G V3, can be regarded as a decorated triangular 
lattice. Its dual is the hexagonal lattice G\ 3 = (V*,El 3 ), where every edge in 
G* 13 is doubled. 

• The subgraph G23 = (V* U V 3 ,E 23 ), where E 23 are the set of edges joining 
vertices e = {ij} with i G V* and j G V 3 , can be regarded as a decorated 
hexagonal lattice. Its dual is the triangular lattice G* 3 = (V,£ , | 3 ), where every 
edge in G 23 is doubled. 

The union-jack and bisected- hexagonal lattices with free boundary conditions 
shown in Figure |3] are quasi triangulations. But in both cases it is very easy to 
modify them slightly so they become Eulerian triangulations: if suffices to add an 
outer vertex belonging to V 3 (black dot) and join this extra vertex to every vertex on 
the outer cycle of each graph. As the degree of all of them was odd, with the new 
edge, it becomes even, and we obtain the desired Eulerian plane triangulation. 

3 Models to be studied 

In this section we will introduce the models we are going to use in the following 
sections, and study their main properties. We cover a well-known model (the AT 
model [H]), a model known in the literature in a different context (the HMBW 
model [2j[12j[T3]), a model used as a mere technical step (but we think it has some 
interest in its own right: the mixed AT model [3]), and a new model (the infinite- 
coupling limit AT model). 

3.1 HMBW model on an Eulerian plane triangulation 

Given an Eulerian plane triangulation G = (V, E) with vertex set V, edge set 
E, and the set F(G) of all the triangular faces, we define the Hintermann-Merlini- 
Baxter-Wu model (HMBW) on G as follows: on each vertex % G V, we place an Ising 
spin <7j = ±1, and these spins interact via the Hamiltonian: 

%HMBW = — J 0-iO~jO~k, (3.1) 

{i,j,k}=teF(G) 
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where the sum is over all the triangular faces t = {i,j,k} of G bounded by the 
vertices i,j,k, and the corresponding 3-spin interaction has a coupling constant J. 
The partition function is 



J HMBW 



(G;J) = ^V Whmbw . (3.2) 
M 



In the ferromagnetic regime J > 0, if G is a plane Eulerian triangulation the 
HMBW model has exactly four ground states: one in which all spins take the same 
value a = +1, and three states in which one spin takes the value a — +1, and the 
other two take the opposite value. 

One important property of the ferromagnetic HMBW model on an Eulerian plane 
triangulation is that this model is self-dual [13] with dual coupling satisfying the 
identity 

v v* = 2 , (3.3) 
where the temperature-like variable v is defined as 



e 2J - 1 . (3.4) 



This self-dual point satisfies: 



v c = V2 (3.5a) 
; 2Jc = 1 + V2 (3.5b) 



In the BW and HM models this self-dual point (13. 5p constitutes the unique critical 
point of the model, as the exact solution shows [2llU[5]. Notice that the critical 
coupling ( I3.5bl) is the same as the critical coupling for the square-lattice Ising model 



3.2 The Ashkin-Teller model 

Let us consider an arbitrary graph G = (V,E), with vertex set V and edge set 
E. The Ashkin-Teller model (AT) is defined on G as follows: on each vertex i G V 
we place two Ising spins 0"i,Tj = ±1, and these spins interact via the Hamiltonian 
( II. 2p [I4"] . Therefore, we can regard this system as two copies of the graph G, with 
the a spins living on one copy, and the r spins living on the other. The a (resp. r) 
spins interact among themselves via a nearest-neighbor coupling K2 (resp. K' 2 ), and 
both copies interact via a 4-spin coupling K4. The partition function for this system 
is 

Z AT {G;K 2 ,K' 2 ,K 4 ) = ^ e^ AT . (3.6) 

The Boltzmann weight for each edge e = (ij) can be read from ( 11.21) /( )3~6l) . and 
depends only on the products OiO^ TiTj = ±1: 

^, ; (a,a r r,7 ; ) = exp [K 2 (Ji(Jj + K' 2 TiTj + K 4 cr i a j T i T j } . (3.7) 
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Indeed, the weights are the same for all edges (ij) G E. For each edge (ij) 6 E, 
there are four possible spin configurations. Their corresponding weights are denoted 
by Uk with k = 0,1,2,3, and they are given in the second column of Table dJ We use 
the standard definitions (see e.g., Ref. [20J): 

u = (3.8a) 

wi = c^(+l,-l) (3.8b) 

w 2 = (3.8c) 

oo 3 = w y (-l,-l) (3.8d) 

We can also use these weights (13. 8p to define the corresponding AT model: 

Z AT (G; K 2 , K' 2 , K 4 ) = Z AT (G; u , u u u 2 , w 4 ) . (3.9) 

The values of the couplings K 2l K' 2 ,K A in terms of the weights w , tui, o; 2 , o; 3 can 
be obtained from the equations given in the second column of Table [TJ 



e 4K ° 


= Uq UJi U) 2 W 3 


(3.10a) 


e 4X 2 




(3.10b) 




UJ 2 003 


e *K> 


wo u 2 


(3.10c) 






e 4X 4 




(3.10d) 




^2 



where K is an arbitrary constant fixing the zero of energy. In Eq. (11.21) we choose 
Kq = for simplicity; therefore the weights uj^ satisfy: 

w uiw 2 W3 = 1. (3-11) 

For our purposes, it is also interesting to compute the normalized weights obtained 
by making the weight corresponding to the configuration cr^ = aj, 7$ = Tj equal to 
one (i.e., Uq = 1). These normalized weights are also displayed in the third column 
of TableU 

The AT model (11. 2p / (13. 6p contains two important particular cases: when K4 = 0, 
we obtain two decoupled Ising models with couplings K 2 and K' 2 . At the other 
extreme, the limit K4 — > +00 corresponds to a single Ising model (er = r) with 
coupling K 2 + K 2 . Finally, the line K 2 = K' 2 = K4, corresponds to the 4-state Potts 
model with Jp tts = ^K 2 . 

The Hamiltonian (II. 2p is invariant under any permutation 7r of the coupling con- 
stants (K 2 , K' 2 , K4), as the fields a, r, <tt play symmetric roles in the model: 

(K 2 , K' 2 , Ki) -> ir(K 2 , K' 2 , Ki), neS 3 . (3.12) 

If the graph G is bipartite, then the vertex set can be split into two disjoints sets 
V = V\ U V 2 , such that each edge e = (ij) satisfies that « 6 VJ, and j = Vj with i ^ j. 
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Configuration 



Weight (Normalized) Weight 





= +1 


' T i T j 


= +1 


U 


_ e K +K2+K' 2 +K 4 




1 


OiOj 


= +1 


> T i T j 


= -1 


Ul 


_ e K +K 2 -K' 2 -K A 


e~ 


-2(K 2 +K 4 ) 


OiOj 


= -1 


) T i T j 


= +1 


UJ 2 


= e K -K 2 +K 2 -K 4 


e 


-2(K 2 +K A ) 


(JiCTj 


= -1 


i TiTj 


= -1 


u 3 


= e K -K 2 -K! 2 +K 4 


e 


-2(K 2 +K! 2 ) 



Table 1: Boltzmann weights for the AT model (II. 2p /f l3T6|) . For each edge (ij) e E, 
we first give the configuration of the corresponding a, r variables, we then quote the 
weight read off from (13. 7p [" Weight"], and the normalized weight obtained by making 
the first one equal to 1 ["(Normalized) Weight"]. K is an arbitrary constant fixing 
the zero of energy for the contribution of each edge (ij) 6 E. In (II. 2p we choose 
Kq = for simplicity. 



In this case, there are additional symmetries because we can flip a, r, or both, on 
any of two vertex subsets V^. Then, the uniform AT model on a bipartite graph G is 
invariant under the transformations: 

(K 2 ,K' 2 ,K 4 ) (-K 2 ,K' 2 ,-K A ) (3.13a) 
(K 2 ,K' 2 ,K A ) (K 2 ,-K' 2 ,-K 4 ) (3.13b) 
(K 2 ,K' 2 ,K 4 ) (-K 2 ,-K' 2 ,K 4 ) (3.13c) 

When K 2 = K' 2 we obtain the symmetric Ashkin-Teller model: 

"HsAT ^ [K 2 ((JiO-j + TiTj) + K-OjOjT,; , . (3.14) 

3.3 The Infinite- Coupling Limit Ashkin-Teller model 

As far as we can tell this model is new in the literature|§ Its motivation is rather 
simple: if we look at the normalized weights in Table HJ we see that in the limit 
K 2 , K' 2 , —K 4 — > +00, with K 2 + K 4 = L 2 and K 2 + K 4 = L' 2 kept finite, one 
has that the normalized weight for o~iOj = TiTj = — 1 tends to zero, while the other 
configurations keep nonzero normalized weights. These weights correspond to the 
infinite- coupling-limit Ashkin-Teller model (ICLAT). The Hamiltonian of this model 
can be alternatively written as: 

%iclat = ~ ^ [{L 2 — V)OiOj + (L' 2 -L)TiTj + \. a, ay , t j , (3.15) 

(ij)dE 



2 After the completion of this work, we learned that a particular case of this model had been 
previously considered by Ikhlef and Rajabpour [21 . Their findings have been summarized in a 
remark in Section 15.11 We thank the referee for bringing this paper to our attention. 
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where we take the limit L — » — oo, with the couplings L 2 ,L 2 kept finite. The corre- 
sponding partition function is 

Z 1C lat(G;L 2j L' 2 ) = lim Ve"^. (3.16) 

L— >-oo z — • 

The Boltzmann weight for a given edge e = (ij) is given by ( 13 . 1 5f) / ( 13.161) and it reads: 

loipiOj, nTj) = exp [(L 2 - L) a^j + (L' 2 - L) + l.o,o , ;, ; , [ ■ (3.17) 

The values of this weight for the four possible spin configuration are given in Table [2j 
Indeed, if we normalize these weights so that the normalized weight for OiOj = r^Tj = 1 
is equal to one (see the third column in Table [2]), we obtain the same weights as those 
obtained by taking the appropriate infinite-coupling limit in the normalized weights 
for the AT model (displayed in the third column in Table [1]). 



Configuration Weight (Normalized) Weight 



OiOj 


= +1 


' T i T j 


= +1 




— gLo+L2+L' 2 


-L 


1 


OiGj 


= +1 


i TiTj 


= -1 


Wi 


— qLq+L2—L 2 


-L 


e -2L' 2 




= -1 


> T i T j 


= +1 




— qLo—L2+L' 2 


-L 


e -2L 2 


(JiCTj 


= -1 


' T i T j 




U 3 


— g-Lo+3L— 1/2- 


~L' 2 


e 2(2L-L 2 -L> 2 ) _^ q 



Table 2: Boltzmann weights for the ICLAT model (13.15P / (I3.16p . For each edge 
(ij) G E, we first give the configuration of the corresponding a, t variables, we then 
quote the weight read off from ( 13 . 1 T[) [" Weight"], and the normalized weight obtained 
by making the first one equal to 1 ["(Normalized) Weight"] in the limit L — > —oo. 
L is an arbitrary constant fixing the zero of energy for the contribution of each edge 
(ij) G E. In (I3.15P we choose L = for simplicity. 



When L 2 — > +oo, then the a spins are all equal; therefore, the model reduces to 
an Ising model on the r spins with coupling L' 2 . Indeed, when L 2 — > +oo, the r spins 
are all equal, and the model reduces to an Ising model on the a spins with coupling 
L 2 . 

In the ICLAT model ( 13. 151) a and r play a symmetric role, therefore, the model 
is invariant under the transformation 

(L 2 ,L' 2 ) -> (L' 2 ,L 2 ). (3.18) 

If the graph G is bipartite, we can define a new spin = = ±1, so that the 
Hamiltonian ( 13. 15ft reads: 

^iclat = - ^ [{L 2 — L)(Ji(jj + (L' 2 - L)piPj(Ti<Tj + Lpipj] , (3.19) 

(ij)eE 
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when L — > — oo. If we flip the spins pi on one of two vertex subsets Vi, then the model 
is invariant under the transformation (L 2 — L, L 2 — L,L) — > (L 2 — L,L — L 2 , —L). We 
obtain 

"Hiclat = ~ ^ [(La - £)0iOj - (A - L)piPj°i°j ~ Lpipj] , (3.20) 

(*i>e-B 

in the limit L — > — oo. By comparing (13. 15)) / ( 13.20p . we see that, if we redefine 
L — > L + L 2 — > — oo, we obtain the equivalent Hamiltonian 

"HlCLAT = ~ ^ [(L 2 - L 2 - Ir)«7iO-j- + (~L 2 - L)piPj + LpiPj<Ti<Tj\ , (3.21) 

in the limit L — > — oo. Therefore, the ICLAT model on a bipartite graph is invariant 
under the transformations: 

(L 2 ,L' 2 ) -> (L 2 -L' 2 ,-L 2 ) (3.22a) 
(L 2 ,L' 2 ) -)■ {-L 2 ,L' 2 -L 2 ) (3.22b) 

where the second equation comes from using the variables ij,p) instead of (a, p). 
If we take L 2 = L' 2 = in (13. 15)) . then it reduces to 

%ICLAT = - 2J [ _ LcTj(Jj - LTjT,- + LViGjTiTj] , (3.23) 

in the limit L — >■ — oo. If the lattice G is bipartite, then we can flip the r and a 
spins on one of the vertex subsets Vi, so that all the couplings become equal to L. So 
it corresponds to the zero-temperature limit of the antiferromagnetic 4-state Potts 
model. 

3.4 Mixed Ashkin-Teller model 

This model appears in a natural way when one consider the duality transformation 
of an AT model. Let us now suppose that G = (V, E) is an arbitrary planar graph; 
hence, it has a planar dual pair G* = (V*, E*), Indeed, one can define an AT model 
on both G or G* as in the preceding section. However, we can define a mixed Ashkin- 
Teller model (mAT) with one set of spins living on G, and the other set on G*, and 
a clever coupling of both Ising models. This coupling is based on the following fact 
about dual graphs: the number of edges is the same in both G and its dual G*, and 
there is natural bijection between these two edge sets: we can always draw G and G* 
in the plane in such a way that each edge e = (ij) G E intersects its corresponding 
dual edge e* = (i*j*) G E* exactly once. 

The formal definition of the mixed AT model is as follows. We place at each 
vertex i 6 V (resp. at each dual vertex i* G V*) of the graph G (resp. of the dual 
graph G*) an Ising spin Oi = ±1 (resp. 7> = ±1), and these spins interact through 
a nearest-neighbor coupling K 2 (resp. K 2 ). These two Ising models are coupled via 
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the four-spin interaction K^OiOjT^Tj*^ where e* = G E* is the unique dual edge 

associated to the edge e = (ij) G E. Therefore, the Hamiltonian for this model is: 

"HmAT = ~K 2 (JiVj - K' 2 ^2 r ** r i* ~ K * a i a o T i* T 3* • ( 3 - 24 ) 

(ij)eE {i*j-)dE* (ij)EE 

The partition function for this mixed model is given by: 

Z mAT (G,G*;ir 2 ,i^,ir 4 ) = Yl e - H ^. (3.25) 

Wt} { r i*} 
i&v v e v* 

It is important to note that the spins living on G couple through the coupling K 2 , 
while those living on its dual G*, couple via K 2 . 

The notation in (13.241) can be lightened by noting again that the correspondence 
(ij) — » (i*j*) is bijective. Thus, we can loosely use the labels ij (instead of i*j*) 
for the t spins (but one should remember that in the mixed AT model, the r spins 
live on the dual vertices V*, while the a spins live on V). We can loosely write the 
Hamiltonian for this model as in Eq. (11.21) . 

Finally, we can also represent the partition function of the mixed AT model (13.251) 
in terms of the weights {tUk} with k = 0, 1, 2, 3 [cf. (13. 8p ] as in Eq. (13.91) : 

Z m AT(G,G*;K 2 ,K' 2 ,K 4 ) = Z mAT (G,G*;u ,u 1 ,u 2 ,u 3 ) . (3.26) 

It is important to recall that, if we interchange G -R- G*, then we should interchange 
the couplings K 2 <-> K 2 , or equivalently, the weights U\ f-> uj 2 . 

Z mAT (G,G*;K 2 ,K' 2 ,K 4 ) = Z mAT (G* ,G; K' 2 , K 2 , K 4 ) (3.27a) 

= Z mAT (G*,G;ui ,ui 2 ,u 1 ,uj 3 ) (3.27b) 



Remark. Note that, contrary to the standard AT model, the number of o spins 
(= \y\) is in general different from the number of r spins (= 

If G = (V, E) is an Eulerian plane triangulation, then it is 3-colorable, and there 
is a natural tripartition of the vertex set V = Vi U V 2 U V3, and edge set E = 
E\ 2 U -E13 U E 23 . If we consider the subgraph G\ 2 = (Vi U V 2 ,Ei 2 ), then its dual 
subgraph is given by G\ 2 = (V3,El 2 ), where the dual edge set E\ 2 is built in the 
standard way. Then, for each edge e G E 12 , there corresponds a dual edge e* G E* 2 
such that it crosses e once. Then the Hamiltonian (I3.24p can be (loosely) rewritten 
as 

UmAT = - 7^ + K 2 TiTj + K A (TiO jTiTj] , (3.28) 

{ij)&E 12 

and the corresponding partition function (13.251) as 

Z mAT (G 12 ,G* 12 ;K 2 ,K' 2 ,K 4 ) = }Z Yl e" WmAT . (3.29) 

i&ViUV 2 i* e v 3 
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We can also take the same infinite-coupling limit as in the previous section, namely 
K 2l K' 2l — K 4 — > +00, while L 2 = K 2 + K 4 and L' 2 = K' 2 + K4 are kept finite. Using 
the above-discussed bijection between E and E*, this model can be written loosely as 
in Eq. f)3. 15j) . but again, keeping in mind that the a (resp. r) spins live on V (resp. 
V*). 



4 Mappings among models 

4.1 Mapping between the AT and the mixed AT models 

In this section we will consider a general AT model with couplings K 2 , K 2 , K4 on a 
planar graph G = (V, E), not necessarily Eulerian. We assume that either K 2 > \K^\ 
or K 2 > \Ki\. Without loss of generality, let us assume that K' 2 > \K±\. The relation 
between the AT and mixed AT models arises when trying to derive the duality relation 
for the former model. We will follow the ideas of Ref. [20J; but instead of looking 
for a relation with the 8-vertex model, we will seek the relation with the mixed AT 
model f l^24j) / d3^5|) . (See also 01221123 ■) 

Let us start with an AT model on a graph G = (V, E) [cf. (11. 2p / (13. 6p ]: 

Z AT {G;K 2 ,K' 2 ,K 4 ) = U eWj ' E II ' ^ 

{a} (ij)GE {*} (ij)GE 

where we have split the sum over the spin configurations into two terms: one over the 
a spins, and the other over the t spins (we use the letter t to denote the second set 
of Ising spins in this case). The effective edge- dependent two-spin coupling between 
the t spins is given by 

K 2ij = K^ + KiOiVj. (4.2) 

Notice that the condition K' 2 > \K^\ implies that K 2i j > for all edges (ij) G E. 
We now do a high-temperature expansion using the fact that 



e 



cosh(K 2ij ) 1 + Utj tanh(i^2ij) • (4.3) 



When we expand the product over the edges e G E, we get a sum over spanning 
subgraphs E' C E, where E' is the subset of edges contributing with the factor 
titj ta.nh.(K 2i j) . When we sum over all possible configurations of the t spins, only 
those subgraphs E' with an even number of incident occupied edges on every vertex 
x £ V survive. In summary, the partition sum reads 

Z AT (G;K 2 ,K' 2 ,K 4 ) = 2l y l^ J] [e K2aiaj oosh{K 2ij ) 

{a} (ij)EE 

j2 n tanh (^), (4.4) 



x 



E' C E (ij)eE' 
E' Eulerian 



where K 2i j is given by (14. 2p . 
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We now introduce the new Ising variables T{ = ±1 living on the dual vertex 
set V*. They are assigned such that for any Eulerian spanning subgraph (V,E') of 

g n •./•;): 

e=(ij)eE' ^ r t ,r r = -1, (4.5) 

where are the vertices of the dual edge e* corresponding to the edge e. (See 

Figure SJ) In words, given an edge e G E, if this edge belongs (resp. does not belong) 
to the subset E' C E, then the r variables associated to this edge have different (resp. 
equal) signs. 




(a) (b) (c) 

Figure 4: (a) Dual edge e* = G E* corresponding to an edge e = (ij) G E. 

The vertices i, j belong to V, and the dual ones i*, j* to V*. (b) When e belongs to 
the spanning subgraph (V,E') (depicted as a thick line), then 7>Tj» = —1 (depicted 
as dots of different colors), (c) When e does not belong to the spanning subgraph 
(V, E'), then t^tj* = 1 (depicted as dots of equal color). 




Figure 5: Mapping (14. 5 p between Eulerian spanning subgraphs on a vertex of degree 
3. The solid black (resp. dashed gray) edges belong to E (resp. E*), and the dots 
correspond to vertices of V*. If e G E belongs (resp. does not belong) to the spanning 
subgraph (V, E 1 ), the corresponding edge is depicted as a thick (resp. thin) solid black 
line. In panels (a)-(d) we show the four different configurations: no edge in E' (i.e., 
all Ising spins have the same value), and two edges in E' (i.e., two Ising spins take 
the same value, and the other one takes the other value). 



For instance, the relation between edge configurations and r spin configurations 
for vertices of degree 3 and 4 are given, respectively in Figures [5] and [6j For a vertex 
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Figure 6: Mapping (I4.5|) between Eulerian spanning subgraphs i on a vertex of degree 
4. In panels (a)-(g) we show the eight different configurations: no edge in E' (i.e., 
all Ising spins have the same value), two edges in E' (i.e., two consecutive Ising spins 
take the same value, and the other two take the other value; or one spin takes one 
value, and the other three spins take the opposite value), and four edges in E' (the 
spin values alternate as we move around). 



v G V of degree d, the number of possible Eulerian configurations (i.e., an even 
number of edges incident to v) is given by 



n>0 



d 
2n 



td-l 



(4.6) 



which is exactly the same number of r-spin configurations of the d neighboring ver- 
tices (modulo a global reversal of the spin values). Indeed, there is a one-to-two 
correspondence between edges subsets E' and {r} configurations: for each subset E' 
there are two equivalent r configurations differing by a global change of sign. 

Note that the a (resp. r) spins live on the vertices of V (resp. V*). Thus, we 
arrive at a mixed AT model defined on G and G*\ 



Z AT (G;K 2l K ; 2: K 4 ) 



2 \y\ 



-i 



E E II [e K ^ cosh(K 2lJ ] 



{n} 



1 + 7>T- 



>< n 

(ij)dE L 

2 \V\-l II UijfaVj'Ti'Tr) 

{vi} {n} {ij)£E 

iev i g v* 



(4.7a) 
(4.7b) 



where K 2 ij is given by ( 14. 2p . and to each edge e = (ij) G E there corresponds a 
unique dual edge e* = G E*, and vice versa. The extra factor 2 _1 comes 
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from the above-mentioned two-to-one relation between r configurations and Eulerian 
subgraphs E' C E. The second term in the product is just a way to express the fact 
that when e £ E' (so -7> = —Tj*), there is a contribution tanh(if2jj); an d when e i?', 
(so Tj* = Tj*), then the contribution is 1. The Boltzmann weights Qij associated to 
the above representation (14.71) are given in Table [3j 



Configuration Weight (Normalized) Weight 

a i a i = +1 , TiTj 

°i°3 = +1 J T i T j 
= ~ 1 > T i T j 
°i°3 = _1 J T i T j 



Table 3: Boltzmann weights for the model (14. 7p . For each configuration of the a, r 
variables, we first quote the weight (up to an unimportant global factor) read off from 
(14. 7p [" Weight"], and then the normalized weight obtained by making the first one 
equal to 1 ["(Normalized) Weight"]. 



+1 uj = e K2 cosh(i^ + K*) 1 

-1 ui = e K2 sinh(K 2 + K 4 ) tanh(iT 2 + K 4 ) 

+1 Q 2 = e-^cosHK' 2 -K 4 ) e---S|pg 

-1 S 3 = e-^sinh(^-K 4 ) e-^^feS 



The relation between the weights Qk of the transformed model (14.71) and the 



original ones Uk is given for every edge 


(ij) G £ by 




U = 


-(uo + uj 


(4.8a) 


U X = 


i(w -wi) 


(4.8b) 


u 2 = 


-(u 2 + u 3 ) 


(4.8c) 


u 3 = 


~(u 2 -us) 


(4.8d) 



The above results can be summarized in the following 

Theorem 4.1 Let G = (V,E) be a planar graph. Then the Ashkin-Teller model 
( II . S\) /( EZ5P on G, with K' 2 > \K^, is equivalent to a mixed Ashkin-Teller model 
tfT£fl )/ (ZM) on G and G* such that 

Z AT (G;uo,ui,u 2 ,u 3 ) = 2 |v1 ~ 1 Z mAT (G, G*\ Q ,ux,u 2 ,u 3 ) , (4.9) 

where the weights u^ are given in terms of the weights u^ by fl^.#| ). 

Remark. Let us remind that in the mixed AT model in Theorem 14. II the a (resp. r) 
spins live on V (resp. V*). The a (resp. r) spins have a nearest-neighbor coupling 
K 2 (resp. K' 2 ); and both a, r spins interact via a 4-spin coupling K4. 
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The duality transformation for the AT model can be easily obtained from Theo- 
rem UJ] by playing the same game. We want to express the partition function of the 
original AT model on G in terms of the partition function of another AT model on 
the dual graph G* . This is easily done as follows: 

2^ y|+1 Z AT (G;w ,Wi,W2,W3) = Z mAT {G,G*;u} , Wi,w 2 ,a5 3 ) (4.10a) 

= Z mAT (G*,G;u} ,u} 2 ,u} 1 ,u} 3 ) (4.10b) 

= 2-W* +1 Z AT (G* ] cu*,u* 1 ,u* 2 ,u* 3 ) (4.10c) 

where the first equality is just Theorem 14.11 the second line comes from using the 
symmetry relation (13.271) for the mixed AT model, and the last one comes from using 
Theorem 14.11 on the dual graph G*. The expression of the new weights ul in terms 
of the original ones is 



* 


= UJ + UJ 2 = 


^(w + u 1 +u 2 + u 3 ) 


(4.11a) 


* 


= uj -uj 2 = 


!/ \ 
-(uio + u 1 -u 2 - u 3 ) 


(4.11b) 


* 

u 2 


= udi + Q 3 = 


!/ \ 

~{U) + UJ 2 -Ui- u 3 ) 


(4.11c) 


* 

u 3 


= Q 1 - Q 3 = 


l , \ 

-{Uq + U 3 -UJi- U 2 ) 


(4.11d) 



The above results can be summarized in the following 

Theorem 4.2 (Wu [20]) Let G = (V,E) be a planar graph. Then the Ashkin- 
Teller model lll.ty) / (37b}) on G is equivalent to another Ashkin-Teller model on the 
dual graph G* such that 

2"l y l Z AT (G; coo, Wi,W2,w 3 ) = Z AT (G*; u*, wj, u* 2 , u* 3 ) , (4.12) 

where the weights ul are given in terms of the weights u>k by ( TjTTT ). 

If G is self-dual, then the corresponding AT model is self-dual when the following 
condition holds: 

U = LUi+LU 2 + U 3 . (4.13) 

This situation occurs for instance for the infinite square lattice. 

4.2 Mapping between the mixed AT and the ICLAT models 

Let us now explore a little further the relation between the mixed AT on G and 
G\ and the AT model on G*. We know from (Onjl / ljOIj l that: 

Z mAT (G, G*;u;o,wi,W2,W3) = 2" |v " l+1 Z AT (G*\ ujq, lu 2 , wjj) (4.14) 
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where the new weights are given by 

Uq = u + u 2 (4.15a) 

u\ = ojq — 0J2 (4.15b) 

u>2 = oj\+ u 3 (4.15c) 

W3 = oj\ — U3 (4.15d) 

It is easy to see that = if and only if u\ = U3. And the last equality holds true 
if and only if K 2 = K4. In this particular case, we have 

2l y *l- 1 Z mAT (G,G*;u; ,u; 1 ,u; 2 ,u; 1 ) = Z AT (G*; < u{, < 0) (4.16a) 



Z AT (G*; 1,^,^,0) 11 ^0 (4.16b) 



u ' fe'>6£; 

But the ratios of these weights are precisely the couplings of the ICLAT model: 

U * „-2KZ'-2KZ . „-2L% 



e~ 2K * - 2K * = e- 2L * (4.17a) 



^ = e -2KS-2Kl = e -2L* 2 (417b) 

In terms of the original couplings, we get that the couplings of the ICLAT model are 
given by the expressions: 

e _ 2L ,' = ^0^2 = tanh(2 ^ 2) (41ga) 

e-^ = = — tt~77T (4.18b) 

Wo + w 2 cosh(2A 2 ) 

The above considerations prove the following 

Theorem 4.3 Let G = (V, E) be a planar graph. Then the mixed Ashkin-Teller 
model \3.2J$ / \3. 25\) on G and G* with equal couplings K 2 = K4, is equivalent to the 
infinite-coupling-limit Ashkin-Teller model on G* , such that 



ZmAT{G, G*] LOq, UJ\, (x>2, (jJ\) — 2 l y l + 1 



n 

(ij)eE 



u 



Z 1CLAT (G*;L;,L^) (4.19) 



where the weights couplings L\,L\ are given in terms of the couplings K^K^ by 



If we prefer to write the partition function of the mixed AT model as an ICLAT 
model on G, the procedure is similar: in this case we have to use (14. 8[) to define the 
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couplings of the AT model ojk in terms of those of the mixed AT model Qk- We obtain 
(compare to (I4.15P ) 

u Q = u + cDi (4.20a) 

uji = Q - ui (4.20b) 

uo 2 = u 2 + o) 3 (4.20c) 

u 3 = u 2 — o) 3 (4.20d) 

so that 

Z mAT (G,G*; uj , u)i, 0)2,0)3) = Z at {G,Uq,ojx,u)2,ojz) (4.21) 

We then notice that 073 = if and only if 0)2 = 0)3, which is in turn equivalent to 
K' 2 = K 4 . 

But the ratios of weights for this AT mode are precisely the couplings of the 
ICLAT model: 



^ = e -2K^-2K t = e -2L> (422a) 
e -2K 2 -2 Ki = e -2L 2 (422b) 



0J 

In terms of the mixed AT couplings we get that the couplings of the ICLAT model 
are given by the expressions: 

e -2i/ 2 = = tanh(2ir 2 ) (4.23a) 

e- 2i2 = J^- = ~ (4.23b) 

wo + wi cosh(2iT 2 ) 

The above considerations prove the following 

Theorem 4.4 Let G = (V, E) be a planar graph. Then the mixed Ashkin-Teller 
model \3.2J^ )/ \3. 25\) on G and G* with equal couplings K 2 = K 4 , is equivalent to the 



infinite-coupling-limit Ashkin-Teller model on G, such that 



Z mAT (G,G*;Lo ,Ld 1 ,Q 2 ,^2) = 2 



n w ° 



Z lCLAT (G;L 2 ,L' 2 ) (4.24) 



where the weights couplings L 2 , L' 2 are given in terms of the couplings K 2 , K 2 by 



4.3 Partial trace transformation 

The goal of this section is to show that the HMBW model on various Eulerian 
plane triangulations can be mapped onto ICLAT models defined on the most common 
lattices (square, hexagonal, triangular, and kagome). Let us start with the following 
lemma: 
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Lemma 4.5 Let a±, . . . , o~2k take the values ±1, and define o~2k+i = o~\. Then, if 
l<k<3, 





2k 




cosh 




= cosh 









k k 
J^<?2i + ■ , y2^2,j 1^2, ^2, 



+ 1 



i=l 



(4.25) 



Moreover, when k > 4 i/iis equality does not hold. 



PROOF. The statement (I4.25P is equivalent to prove that the following quantity A(k) 
vanishes when all of = 1, o 2 fc+i = o"i, and 1 < k < 3: 



A(A;) 



2A: 



2/c 



,i=i 



5^0"2i + ^ 0-2j-lO- 2 jO-2j+l 



(4.26) 



The most direct way to show that A(k) = for k = 1,2,3 when of = 1 and 
d"2fc+i = o~i is to write down the above expression with U = 1 — of and taking into 
account the latter condition (o 2 £ + i = 01). We arrive easily at the expressions: 



A{1) = tl(t 2 -l), 
A(2) = -tit 3 (o- 2 + o 4 ) 2 



(4.27a) 
(4.27b) 



Aft) = - ^^hj-i [hj+ihj - 2o 2 j0 2 j_2 + ^j^j +1^-2^-3] , (4.27c) 
3=1 

where we have explicitly used the identifications o 2 fc+i = 01 and 02k = o~o- It is 
obvious that when ti = 0, then A(k) = for k = 1,2, 3. 

For k > 4 ( I4.26P cannot be true as the first term of the r.h.s. of Eq. (I4.26P 
contains products of up to four spins, while the second term of the r.h.s. of Eq. (I4.26P 
contains products of six spins, and none of them can be simplified using the boundary 
conditions o 2 fc+i = o~\- ■ 



Remark. Eq. (I4.25P can be rewritten as a sum over a spin o : 

2k 



Y e i=1 = ex p 



<T0=±1 



(T =±l 



k k 

j^o o 2i + j yVoP 2 j~io 2 jP 2 j + i 



1=1 



. (4.28) 



Therefore, a 3-spin Ising model can be written as a certain Ising model with 2- and 
4-spin interactions. 

The first non-trivial application corresponds to k = 2 (see Figure [7J). Let us 
suppose we have an Eulerian plane triangulation G = (V,E), such that V = V\ U 
V 2 U V3, and such that each vertex i e V 3 has degree A 3 = 4. (This happens for 
all the lattices belonging to the family described in Section 12.21 ) If we perform 
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(a) (b) 

Figure 7: (a) The Ising model with uniform 3-spin interaction J{aoa\a 2 + o-Q(j 2 a 3 + 
aoa 3 a^ + aoa^ai). (b) The Ising model with 2-spin interaction J{aoa 2 + cr a4), and 4- 
spin interaction J(ao<Jia 2 a 3 + <7o0"30"40"i) obtained from (a) after applying Lemma I4T51 
with k = 2. 

the transformation of Lemma 14.51 to the vertices on V3, then the HMBW model on 
G with coupling constant J (13. ip / (13. 2p can be written as a mixed AT model on 
G13 = U V3, Si 3 ), whose dual is G* 13 = (V 2 , E* 3 ) [or on G 23 = (V 2 U V3, £ 23 ) whose 
dual is G 23 = (Vi,E% 3 )]. Both the two-spin coupling <7i<7j, and the 4-spin coupling 
OiOjO{*Oj* (for any (ij) G -E13) take the same value J; but the coupling between dual 
spins is zero. Again, the edge e* = G E{ 3 is the one uniquely associated to 

the edge e G £13. Notice that the vertices in V 3 (depicted as black dots in Figure [7]) 
subdivide the edge joining the two vertices in V\ (depicted as gray dots in Figure |7J). 
Therefore, there are two dual edges (depicted as blue curves on Figure [3(b)) joining 
the two dual vertices in V 2 (depicted as white dots in Figure [7]). The Hamiltonian of 
the transformed system is therefore 

% = —J (TiO-j — J OiOjOi*Oj* (4.29a) 

= -J ^2 ( X + o>o"j*) (TiOj (4.29b) 

The relation between partition functions is clear: given an Eulerian plane triangula- 
tion G = (G, V), such that V = V x U V 2 U V 3 , E = E 12 U E 13 U E 23 , and the degree of 
all vertices in V3 is A3 = 4, then 

Zumbw(G', J) = Z m Ar(Gi 3 , G* 13 ; J, 0, J) , (4.30) 

where G± 3 = (Vx U V3, Ei 3 ), and its dual graph is G{ 3 = (V 2 , E\* 3 ). 

The k = 3 case is shown in Figure El Let us suppose we have an Eulerian plane 
triangulation G = (V,E), such that V = V\ U V 2 U V 3 , and such that each vertex 
% G V 3 has degree A3 = 6. (This happens for the triangular and bisected- hexagonal 
lattices.) If we perform the transformation of Lemma [4.51 to the vertices on V 3 , then 
the HMBW model on G with coupling constant J (13. ip / (13. 2 p can be written as a 
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0~2 0"5 




Figure 8: (a) The Ising model with uniform 3-spin interaction J{aoa\a 2 + (Togo's + 
0W3O4 + cxoCT4(T5 + 0"oO"5<T6 + a^a^ai). (b) The Ising model with 2-spin interaction 
J(o~oa 2 + croO"4 + °~oO~6), and 4-spin interaction J(aoaia 2 a 3 + OQa 3 a^o^ + 0"o0"50Wi) 
obtained from (a) after applying Lemma [4.51 with k = 3. 



mixed AT model on G% 3 = (Vi U V 3 ,Ei 3 ) [or on G 23 = (V 2 U V^E^)], whose dual 
graph is = (V 2 , •£'13); and again all couplings are equal to J, except the couplings 
between dual vertices, which are zero. The model can be written as in (j4.29p /( |4.30p . 
We can gather the above results in the following 



Theorem 4.6 Let G = (V, E) be an Eulerian plane triangulation with vertex set 
V = V\ U V 2 U V 3 , edge set E = E 12 U E Vi U E 2 z), and such that the degree of all 
vertices in V 3 is A3 = 4 or 6. Then the Hintermann-Merlini-Baxter-Wu model 
on G Ii3. l\) / l[3~2]) with coupling constant J is equivalent to a mixed Ashkin-Teller 
model {3.2J$ / \3.2ty with uniform coupling constant K 2 = K4 = J and K 2 = on 
G 13 = (Vi U I/3, E 13 ) [with G* 13 = (V 2 , E{ 3 )], such that 



■"HMBW 

(G; J) — -^mAT(Gi3, G\ 3 ] J, 0, J) . 



(4.31) 



(We can equally define the mixed AT model on G 23 = (V 2 U V 3 , E 23 ), and the formulas 
are the same after interchanging 1 •<->■ 2 .) 

Using Theorem 14. 3[ we can obtain the following 

Corollary 4.7 Let G = (V, E) be an Eulerian plane triangulation with vertex set 
V = V1UV2UV3, edge set E = Ei 2 U Ei 3 U E 23 ) , and such that the degree of all vertices 
in V3 is A3 = 4 or Q. Then the uniform Hintermann-Merlini-Baxter-Wu model on 
G Ii3. l\) / nT2)) with coupling constant J is equivalent to an infinite-coupling-limit 
Ashkin-Teller model HlBf / HlM on G* n = {V 2 ,E* n ), where G 13 = {V x U V 3 ,E 13 ), 
and such that 



J HMBW 



(G;J) = 2-M +1 Z lCLAT (G* 13l L 2 ,L' 2 ), 



(4.32) 
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where 



e ~ 2L2 = uo n ( 433a ) 
cosn(2J) 



e" 2L 2 = tanh(2J) (4.33b) 



These couplings satisfy 



e -2L 2 \\{ e -2L>\ 2 = l (434) 



We can apply this corollary to the main lattices in this paper: 

(a) If G is the triangular lattice and we choose G13 to be one of its hexagonal sub- 
lattices, then G[ 3 is the other triangular sublattice. Then, the triangular-lattice 
HMBW model with coupling J is equivalent to a triangular-lattice ICLAT model 
with couplings (L 2 , L' 2 ) satisfying the condition (14 .34 p . 

(b) If G is the union-jack lattice and we choose G13 to be a decorated square lattice 
(formed by the degree-4 black vertices and the degree-8 gray vertices in Fig- 
ure E^a)), then G* l3 is also a square lattice, tilted 45°, formed by the degree-8 
white dots in Figure |^a), and with double edges. Then, the union-jack-lattice 
HMBW model with coupling J is equivalent to a square-lattice ICLAT model 
with couplings (L 2 , L' 2 ) satisfying condition 

e -2L> 2 + e -2L' 2 = ^ ( 4 35) 

because of the double edges. 
If G is the bisected-hexagonal lattice, then we have several choices for G13: 

(c) If G 13 is the decorated triangular lattice formed by the degree-4 and degree-12 
vertices in G (black and gray dots, respectively in Figure [3](b)), then G* 3 is a 
hexagonal lattice formed by the degree-6 vertices in G (white dots) and dou- 
ble edges. Then, the bisected-hexagonal-lattice HMBW model with coupling 
J is equivalent to a hexagonal-lattice ICLAT model with couplings (L 2 ,L' 2 ) 
satisfying the condition (I4.35p . 

(d) If G13 is the hexagonal lattice formed by the degree-4 and degree-6 vertices in 
G (black and white dots, respectively in Figure [3](b)), then G{ 3 is a triangular 
lattice formed by the degree-12 vertices in G (gray dots), and double edges. 
Then, the bisected-hexagonal-lattice HMBW model with coupling J is equiva- 
lent to a triangular-lattice ICLAT model with couplings (L 2 , L 2 ) satisfying the 
condition (I4.35p . 

(e) If G13 is the diced lattice formed by the degree-6 and degree-12 vertices in 
G (white and gray dots, respectively in Figure Elb)), then G{ 3 is a kagome 
lattice formed by the degree-4 vertices in G (black dots). Then, the bisected- 
hexagonal-lattice HMBW model with coupling J is equivalent to a kagome- 
lattice ICLAT model with couplings (L 2 ,Z/ 2 ) satisfying the condition (I4.34p . 
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In conclusion, the HMBW model on several Eulerian triangulations can be exactly 
mapped along certain curves ( (14.341) or (I4.35P ) on the ICLAT model defined on other 
lattices. 

Some of these results can be obtained using the decimation transformation (see 
Appendix^]). This latter transformation can be applied in the case k = 2 with general 
couplings Jf, but reduce to some of the above results in the homogeneous case J, = J. 



4.4 Mapping between the BW and two-color non-overlapping 
Eulerian-bond models 

Let us start with a BW model: that is, a HMBW model fETO / (ET2)) model defined 
on the triangular lattice G = (V, E). Then, both the vertex and the edge sets can be 
partitioned into three disjoint sets V — V\ U V<i U V3, and E = E\ 2 U £13 U E 2 z, so that 
for all i 7^ j, Eij contains all edges e = (xy) with x G Vi and y G Vj. Each triangular 
face t = {x,y, z} contains one edge from each set E^; therefore, we can rewrite the 
Hamiltonian (13. ip as 

"Hbw = -J 22 (y i <J A (y p + a q) 7 ( 4 - 36 ) 

{i,j}&E 23 

where {i,j, p} and {i,j, q} are only two triangular faces sharing the common edge 
e = (ij) G E23. Notice that the subgraph (V 2 U V 3 , E23) is a hexagonal lattice. 

The Boltzmann weight associated to an arbitrary edge e = (ij) G E 2 s can be 
rewritten as 

e Ja t a 3 {a p +a q ) = ^ rj^ + ^ ^ + ^ tanh rj^ + ^j} _ 

If we now use the trivial identities: 

, r N1 cosh2J lcosh2J if a v = a a 

cosh [J(a p + a, = - - = \ p q (4.38a) 

(cosh 2 Jy a v - a i [1 if a p ^ a q 

, , 1 - , , 1 CTn tanh 2 J if <j„ = cr a , , \ 

tanh [J(a p + a q )] = a p S ap>aq tanh 2 J = 1 p p q (4.38b) 

1 il o v / a,, 

we can perform a high-temperature expansion of the partition function (13. 2p /f l4~36]) : 
Zbw(G;J) = (cosh27)« 2W«I £ £ n (J_Y'*^ 

{<J V } E> C £ 23 ee£; 23 

p £ Vi -E' Eulerian 

x II vW 9 tanh2J, (4.39) 

eg-E' 

with E' being the subset of edges in E23 contributing with a factor proportional to 
tanh 2 J. Notice that 
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(a) The sum over subsets E' only contains Eulerian subgraphs (V2UV3, E'). Because 
(V2 U V3, -E23) is a hexagonal lattice of maximum degree A = 3, the subgraphs 
that contribute to the partition function (I4.39P are precisely the class of non- 
intersecting Eulerian loops. 

(b) The spins that appear in (14.391) live on the sublattice V\. Because of the factor 
c p S ap}(Tq , all spins on any boundary of a loop should be equal. 

(c) Because the hexagonal lattice is bipartite, all loops have even length. Therefore, 
the factors a p in (I4.39P cancel out, and always give a contribution equal to 1. 

With these considerations (I4.39P can be rewritten as 

Z BW (G;J) = (cosh2J) |i?231 £ (tanh2 J) 1 *' 1 

E' C E 23 
E' Eulerian 

p e Vi 

The inverse powers of cosh 2 J can be interpreted as contributions from domain- 
wall loops: i.e., low-temperature loops which separate regions with all spins equal 
to +1 from regions with all spins equal to —1. Notice that (cosh 2J)~ 5<t j"~ <t 9 only 
gives a non-trivial contribution when a p 7^ a q . These low-temperature loops are also 
non-intersecting, and in addition, they cannot cross any high-temperature loop. See 
Figure [9] for an example. 

In summary, the partition function of the BW model (13. ip / (I4.36P is equivalent to a 
loop gas with the following two conditions: (1) the loops are non-intersecting, and (2) 
there are two species of loops: the high-temperature loops with weight tanh2J per 
edge, and the low-temperature loops with weight cosh -1 2 J per edge. The partition 
function can be finally written as 

k ( 

Z BW (G;J) = (2cosh2J) lE ' 23l J2 f[\ (tanh2J) ij + 

k>0 loops with k j =1 I 
components 

(4.41) 

where Ej is the length of the j-th loop. 

Because tanh 2 J + cosh - 2J = 1, we can parametrize the weights with a single 
parameter 6: 

sin 9 = — I , cos# = tanh2J, (4.42) 

cosh 2 J 

so that 

k 

Z BW (G;J) = (2cosh2J) lE2il Y, Yl Tl i sin ^ 9 + cos ^} • ( 4 - 43 ) 

k>0 loops with k i = l 
components 



cosh 2 J 
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Figure 9: The BW model as a two-color non-overlapping loop model on the hexagonal 
lattice. We show a finite subset of a triangular lattice G = (V,E). The hexagonal 
sublattice is formed by the subgraph (V2 U V3, E23)] the vertices of this sublattice are 
not depicted, while the edges are depicted as solid lines. The vertices belonging to the 
subset V\ are depicted as dots: a black (resp. gray) dot corresponds to an associated 
spin taking the value +1 (resp. —1). The edges not in E 2 z are depicted as dashed 
thin lines. The low-temperature loops (depicted as thick black lines) separate regions 
with opposite spin values. The high-temperature loops (depicted as thick solid gray 
lines) are surrounded by spins all taking the same value. As this is a finite piece of 
the triangular lattice, we can only show some parts of the loops. 



If we perform a duality transformation in the HMBW model with J — > J* (v — > 
2/f *) [13] . the loop weights are interchanged: tanh 2 J — > cosh -1 2 J*, and cosh -1 2 J — >■ 
tanh2J*. This agrees with the interpretation of those loops as high-temperature and 
low-temperature loops. The uniform HMBW model is self-dual when J = J*. In 
terms of the loop model, we obtain that (14.431) is self-dual at 

sin6» c = cos# c = -^=. (4.44) 
Notice that the loop model (I4.43P is an inhomogeneous 0(2) model [24J 

k / n \ 

zomWu ...,&) = £ H E 3 " ( 4 - 45 ) 

fc>0 loops with k 3=1 \a=l / 
components 

with weights 

/3i = tanh2J, fa = cosh -1 2J. (4.46) 

This model becomes homogeneous at the self-dual point (14.441) . and this self-dual 
point corresponds to the critical point for the 0(2) model [23] on the hexagonal 
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lattice: 

c (n) = - 1 (4.47) 

as 

lim0 c (n) = -= = tanh2J c . (4.48) 

This is indeed the critical value Baxter and Wu found for the BW model [cf. ( I3.5bl) . 

Remarks. 1. It is worth noticing that when (3i — > oo and (3% — >■ 0, we obtain one 
type of fully packed loops (FPL) [25]. This model is equivalent to a dimer covering 
of the hexagonal lattice (26], which in turn can be described in the continuum limit 
by a conformal field theory of central charge c = 1 J27JI2H]- However, no value of J 
seems to provide these limiting values. 

2. There exists a direct mapping between the ICLAT model and the two-color non- 
overlapping Eulerian-bond model. The loops correspond to low-temperature graphs 
(or domain walls) of the two copies of Ising spins. We can rewrite the ICLAT-model 
partition function (13. 161) as 



+ (l - e- 2L * - e~ 2L >) <WA,J -(4-49) 



{cj,t} e=(ij)EE 



Indeed, the weights for each edge coincide with those normalized weights given in the 
third column of Table |2j One places occupied bonds on the edges of the dual graph 
G* as in the standard low-temperature expansion procedure: place an occupied bond 
of color "A" (resp. color "B") on the dual edge e* if the Ising spins r (resp. a) on the 
corresponding original edge e = (ij) are not equal, i.e. t^t, = —1 (resp. a^aj = — 1). 
Note that any dual edge cannot be occupied by two bonds of distinct colors, as the 
configuration TjTj = a^Oj = — 1 has zero weight. If G is the triangular lattice, then 
G* is the hexagonal lattice, and loops cannot intersect. But for other lattices, loops 
of distinct colors can intersect at certain vertices. In summary, one obtains a one- 
to-two correspondence between a two-color non-overlapping Eulerian-bond model on 
the dual graph G* , and the ICLAT model on the original graph G. 



5 Phase diagrams 

In this section we will qualitatively describe the phase diagrams for the ICLAT 
models defined on the square, triangular, hexagonal, and kagome lattice. We will use 
the analytic results obtained in the preceding section, as well as exact results from 
the 8-vertex model [3]. 

Some important points will be studied using Monte Carlo (MC) simulations. In 
these cases, we have used the embedding algorithm for the AT model [2H] with one 
difference. In the ferromagnetic regime, in the algorithm of Ref. [21] one fixes one type 
of spins (either a or r) and performs a standard Swendsen-Wang cluster algorithm 
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[3U|I3T] with the other type. In the simulation for this paper, we fixed one of the three 
types of spins (a, r, or err). The fixing of the variable err seems important in reducing 
the critical slowing down for the ICLAT model. For the antiferromagnetic regime, we 
performed the embedding using the Wang-Swendsen-Kotecky [32] algorithm. 

The Binder Qi cumulants for different observables were computed for each value of 
the parameters {L%, L' 2 ), and from the crossing of these cumulants for distinct values 
of the linear size, we extracted the critical values of the parameters. Details will be 
published elsewhere. 



5.1 Square- lattice ICLAT model 

In this section we will consider the phase diagram for a square-lattice ICLAT 
model. It is interesting to start with the phase diagram for the symmetric AT model 
on the square lattice, as it is simpler than that of the full AT model [33]. 

First of all, as the square lattice is bipartite, the phase diagram of the symmetric 
AT on the square lattice should be invariant under the interchange K 2 — > —K 2 [cf. 
f)3.13p ]. Thus, we will focus on the part with K 2 > 0: See Figure [TUl 

There are some "easy" points in this phase diagram. 

• The line K4 = corresponds to two decoupled Ising models, so there are Ising 
critical points at (K±, K 2 ) = (0, ±| log(l + v2))- Point DIs in Figure [10] repre- 
sents the ferromagnetic one (+ sign). 

• The line K 2 = corresponds to an Ising model in the variable err. There are 
Ising critical points at (K 4 , K 2 ) = (±|log(l + a/2),0). Point Is in Figure [TUl 
represents the ferromagnetic one (+ sign), and point AFIs represents the anti- 
ferromagnetic one (— sign). 

• The limit — > 00 corresponds to an Ising model with coupling 2K 2 ; there- 
fore, we have Ising critical points at (K4, K 2 ) — > (oo,±|log(l + a/2)). The 
ferromagnetic point (+ sign) is represented by point Is' in Figure [TUl 

• The K 2 = K4 subspace corresponds to the 4-state Potts model. This model 
has a ferromagnetic critical point at [K^K 2 ) = ^log3 (point P in Figure [TUl . 
The antiferromagnetic 4-state Potts model on the square lattice is disordered 
even at zero temperature [34"[I3"5]. 



The square lattice is self-dual; therefore the square-lattice AT model is self-dual 
on the curve (14. 13ft . For the symmetric AT model (13.141) . this curve takes the form: 

e~ 2Ki = smh(2K 2 ). (5.1) 

The AT model on any planar graph can be mapped onto an 8-vertex model on 
the medial graph [20]. In particular, the AT model on the square lattice can be 
mapped onto a staggered 8-vertex model on the square lattice (which has not been 
solved in general). As a special case, the AT model on the self-dual manifold (14.131) 
maps onto a homogeneous eight-vertex model, which is exactly solvable [3j. Finally, 
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Sq— lattice Ashkin— Teller model 
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Figure 10: Phase diagram on the square-lattice symmetric AT model (I3.14p in the 
plane (K^K 2 ). The self-dual curve (15 .ip is B-DIs-P-C. The solid (blue) curves rep- 
resent the critical part of the self-dual curve, the (red) dot-dashed curves represent 
Ising-like transition lines, the dashed (orange) lines represent the 4-state Potts sub- 
space, and the dotted curve represents the noncritical part of the self-dual curve (15. ip . 
The Roman numerals designate the different phases of the model (see text). 



the symmetric AT model (13.141) on the square lattice maps (after a simple further 
transformation) onto a homogeneous 6-vertex model. In this way, Baxter showed 
that the self-dual curve (15. ip is critical only for K± < | log 3 (blue solid curve in 
Figure [TUT) . and it is noncritical for K A > |log3 (dotted curve in Figure [TUj) . The 
critical part can be described by a conformal field theory of central charge c = 1 (i.e., 
it can be related to the Gaussian model). Along this line the critical exponents vary 
continuously, thus violating the usual notion of universality. 

Even though the exact solution of the symmetric square-lattice AT model is not 
known, we have evidence [36T438J that two critical curves emerge at the 4-state Potts 
model critical point P: one goes to the Ising critical point Is, and the other one 
tends as K 4 — > oo to the critical Ising point Is'. Finally, there is another critical 
curve emerging from the Ising critical point AFIs and pointing toward K± — > —00. 
The exact location of these three curves is still an open problem, as well as their 
universality classes. However, it is generally believed that the three curves belong to 
the Ising universality class. 

The phase diagram of this model shows four different phases: 

I. This is the so-called Baxter phase [36]. The spins a and r are independently 
ferromagnetically ordered. There are four extremal infinite- volume Gibbs states, 
one for each choice for the signs of (a) and (r). The sign of (err) is that of (a) (r). 
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II. 



This is the paramagnetic phase, in which the three spins a, r, and err are 
disordered. There is a unique infinite-volume Gibbs state. 



III. In this phase both cr and r are disordered, but err is ferromagnetically ordered. 
There are two extremal infinite-volume Gibbs states, one for each choice of the 
sign of (err). 

IV. This is the antiferromagnetic analogue of III: cr and r are both disordered, but 
err is antiferromagnetically ordered. There are two extremal infinite-volume 
Gibbs states. 

The critical exponents along the critical part of the self-dual curve can be obtained 
by relating the AT model to the 8-vertex model or to the Gaussian model [3"Tll3"5H4"2"] . 
We can parametrize the critical part of the self-dual curve (15. ip by using the parameter 
/i: 

e 4i ^ = l + 2cos/i, 0<fi<— . (5.2) 

3 

This parameter ft is related to the coupling constant g of the Gaussian model by [42^ 

^ = TT (i - 1) , \<g<±- (5-3) 



Then, using renormalization-group arguments [37] 139] one finds that the critical 
exponents along the critical part of the self-dual curve (15 .ip are given in terms of fi 
as 



2-2/ 



a = 3^fy (5 ' 4b) 
P = ~T~~~~~\ (5-4c) 
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(5.4a) 



7 = -r (5.4d) 

1 4(3 -2y) K } 

P' = t (5.4e) 

4(3 -2y) v ' 

i = J~ A \ (5-4f) 



where y is the renormalization-group eigenvalue related to // via 



2u 4 . , 

y = — , 0<y<-, (5.5) 

7T O 



3 Our g is that of Saleur [42] and equals 2tt times of the K of Kadanoff and Brown [39], and 
Yang [41]. 
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and the magnetic exponents 0, 7 (resp. correspond to the a or r (resp. or) 

magnetization. The value of the critical exponents (15. 4p vary along the self-dual 
curve (j5.ip ; but there are two exponents that remain constant on this curve 



V 
6 



2 - 



7 



1 

4 
15 



(5.6a) 
(5.6b) 



Finally, it is worth commenting that the AT domain walls (considered as extended 
curves) have also critical exponents (of the watermelon type). This has recently been 
investigated by Picco and Santachiara [Sl[44j, and by Ikhlef and Rajabpour [2"T] . 
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Figure 11: Phase diagram on the square-lattice g-state Potts model in the plane (q, v), 
where v is related to the coupling constant J by (15.71) . The upper solid (blue) curve 
Vp = +y/q corresponds to the ferromagnetic critical curve. Its continuation (depicted 
as a dashed - red - curve) corresponds to the BK line vbk = —y/Q- The other two 
branches vaf+ = — 2 + y/4 — q and v af- = — 2 — -y/4 — q correspond respectively to 
the critical antiferromagnetic and its dual counterpart curves. The dots • represent 
the critical values of the model for integer values of q. The square ■ represents the 
intersection of the BK curve t>BK = —\/q with the zero-temperature limit in the 
antiferromagnetic regime v = —1: i.e. the point q*(sq) = 1. 



There is a nice relation between the symmetric square-lattice AT model on the self- 
dual curve (15. ip and the ferromagnetic g-state Potts model at its critical temperature 
[4"5] . Figure [TT1 shows the phase diagram of the g-state Potts model in the plane (g, v), 
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where v is given in terms of the Potts model coupling constant Jp tts b)0 

v = e Jpotts - 1 . (5.7) 

There are four curves where the exact free energy is known [3|H6] : 

v = ±y/q (5.8a) 
v = -2 ± y/A-q (5.8b) 

The curve ()5.8aL ) Vp(q) = yfq is known to correspond to the ferromagnetic critical 
point for this model, and it is self-dual. Its analytic continuation into the Berker- 
Kadanoff (BK) phase pT?T - H9] corresponds to the BK curve (|5.8aL ) ubk(?) — — y/q, 
and it is also self-dual. Curve (!5.8bL ) vaf+(q) = —2 + \/A — q is expected to give 
the critical curve for the antiferromagnetic model; in particular, it gives the exactly 
known values of the critical temperature for q = 2,3. This curve is dual to the 
curve (15.8bL ) vaf-(q) = — 2 — y/4 — q. The curves vaf± bound the BK phase. The 
renormalization group flow is repulsive close to this boundary, and it is attracted by 
the curve v^k lying hi between. This curve intersects the line v = —1 (corresponding 
to the zero-temperature Potts model in the antiferromagnetic regime) at a point 
<?*(sq) = 1. 

The above mentioned relation between these two models stems from the fact that 
both models can be written as six- vertex models j3] . Thus, we can use q to parametrize 
the critical part of the self-dual curve (15. ip . Indeed, this q should be interpreted as 
an effective number of states]^ The result is 

y/q = 2 cos /i, qe [0,4], (5.9) 

where \i is related to the AT couplings via (15. ip / (15. 2p . 

The critical 4-state Potts model corresponds to \i = 0, y = 0, g = 4, and yfq = 2 
(point P in Figure [10]); the two decoupled Ising models correspond to the point 
fx = 7r/2, y = 1, g = 2, and y/q = (point DIs in Figure [10]). The part of the AT 
self-dual curve connecting these two models (y < 1) is characterized by a > 0, and 
corresponds to going from q = 4 down to q = along the ferromagnetic critical curve 
vp. Therefore, < \i < n/2 corresponds to the Potts-model critical curve vf. 

If we keep moving to higher values of \x (i.e., 27r/3 < fi < 7r/2), then we enter the 
BK branch v-qk- This corresponds to y > 1, and yfq < 0. In particular, the specific 
heat is not divergent in this regime (i.e., a < 0), as the energy is not a relevant 
operator in the BK phase. We can move down the curve vbk until we enter the 
unphysical regime of the Potts model: i.e., v = —1, which corresponds to the zero- 
temperature limit of the antiferromagnetic Potts model. We expect that this point 

4 We have used the same letter as for the temperature-like variable in the HMBW model (|3.4|) . 
It should be clear from the context which v we are referring to. 

5 Please note that this not imply that the AT and Potts models are equivalent. The AT model 
is untwisted, so that (unlike the Potts model) the central charge is c = 1 throughout the critical 
regime. 



36 



g* is lattice-dependent, as universality usually does not hold in the antiferromagnetic 
regime. This limiting point along the self-dual curve (15. ip as K 4 — > — oo (point B in 
Figure [TUT) corresponds to \x = 27r/3, y = 4/3, g = 4/3, and y/q = — 1. Its universality 
class is characterized by the following critical exponents (15 .4p : 

u = 2, a = -2, P = \, 0' = \, l=l> V=|, (5-10) 
and the two critical exponents r\ and 5 as in (15. 6p . 



Sq-lattice ICLAT model 




exp(-2L 2 ) 

Figure 12: Phase diagram on the square-lattice ICLAT model (I3.15P in the plane 
(e~ 2L2 , e _2L2 ). The (red) solid curves are Ising-like critical curves, except at the 
dots of coordinates (1/2,1/2), (1,2), and (2,1). Point B of coordinates (1/2,1/2) 
corresponds to the point where the self-dual curve of the symmetric AT model (15.11) 
hits this plane, and the other two points are equivalent to the former one due to 
the symmetries of this model. All these points are depicted as solid circles (•). 
Point A (A) shows where the branch emerging from point AFIs in Figure [10] hits the 
ICLAT model. The dashed curve going through point B corresponds to the union- 
jack HMWB model subspace (14.351) . The squares (□) on the axes correspond to the 
exact Ising critical points located at v2 ± 1. The distinct phases are denoted by 
the ordering of the distinct spins: e.g., (cr)p, (t)f, and (<jt)af means respectively, 
that the a spins are disordered, the r spins are ferromagnetically ordered, and that 
the product err is antiferromagnetically ordered. F corresponds to the ferromagnetic 
phase: (u)f, ( t )f, {& t )f', and P corresponds to the paramagnetic phase: (cr)p, (t)p, 
(or) P . 

Let us now consider the square-lattice ICLAT model. Its phase diagram is shown 
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in Figure [121 an d it is symmetric under the interchange L2 -H- L 2 because of ( 13 . 1 8[) . 
Therefore, we will focus on the part with L 2 > L 2 . 

First of all, we want to know where the self-dual curve ( 15. ip hits the plane 
(e _2Z/2 , e~ 2L 2) in the limit K4 — > —00. We can rewrite ( 15. ip as 

in the limit K4 — > —00 and K% = K' 2 — > +00, with K2 + K4 = L2 and K 2 + K4 = 
L' 2 = L 2 held fixed. Thus, the self-dual curve (15. ip hits the ICLAT model at the 
symmetric point 

i^- 2Li ) - G>£) ■ <» 2 > 

This point is depicted as B in Figure IT2"| and it is equivalent to point B' of coordinates 

(e- 2i2 ,e- 2i2 ) = (2,1), (5.13) 
because of the symmetry (I3.22p . 

Remark. Ikhlef and Rajabpour [211 Section 4] have consider the square-lattice 
ICLAT model at the symmetric point ( 15. lip / (I5.12p . They found that at this par- 
ticular point, this ICLAT model can be mapped onto an integrable 19-vertex model. 
Furthermore, they also mapped this particular point of the square-lattice ICLAT 
model onto an integrable dilute Brauer model. 

In the limit L' 2 — > 00, we have a square-lattice Ising model on the a spins. There- 
fore, we have the following two critical points: 

( -2L 2 e -2L' 2 \ = ± 1,0) , (5.14) 

where the '+' (resp. '— ') sign corresponds to the antiferromagnetic (resp. ferromag- 
netic) Ising critical point. Using ( I3.22p . we see that these two points are mapped to 
points at infinity with slopes a/2 =1= 1. All these points belong to the Ising universality 
class. 

The point L 2 = L' 2 = corresponds to a zero-temperature 4-state Potts antifer- 
romagnet. This model is disordered on the square lattice. Therefore, the point 

(e- 2L2 ,e- 2L ' 2 ) = (1,1) (5.15) 

belongs to the paramagnetic phase, which is indicated with the sign P in Figure [12J 
The point L2 = L 2 — > +00 corresponds to a zero-temperature point in the ferro- 
magnetic regime; therefore, the region around the origin in the (e _2L2 , e~ 2Ij2 ) plane 
corresponds to a ferromagnetically ordered phase, denoted by F in Figure [£2J 

We know from Section 14.31 that the union-jack-lattice HMBW model with cou- 
pling J is equivalent to an ICLAT model on the square lattice on the line (I4.35p . This 
line is depicted as a (black) dashed line in Figure [121 Indeed, the self-dual point for 
this HMBW model corresponds to the critical point (15. 12[) . 



38 



Baxter has shown that the AT model on the square lattice can be mapped onto 
a staggered eight- vertex model on the square lattice [3], with weights (a,b,c,d) on 
one sublattice, and weights (a, b, d, c) on the other sublattice. These weights can 
be written (up to some unimportant multiplicative constant) in terms of the pa- 
rameters L 2 , L' 2 of the square-lattice ICLAT model as follows (see [31 Eqs. (12.9.6) 
and (12.9.17)]): 

a = Uq + Ui = l + e~ 2L ' 2 (5.16a) 
b = LU2-LU3 = e~ 2L2 (5.16b) 
c = u 2 + uj 3 = e" 2La (5.16c) 
d = wo-wi = l-e- 2L 'i (5.16d) 

where the weights Uk correspond to the ICLAT model, and can be read from the third 
column of Table HJ The above staggered eight-vertex model becomes homogeneous 
when c = d. In terms of the couplings L,2,L 2 , this implies that the model should 
satisfy (|4.35p . The vertex weights (15 . 1 6[) reduce along this line to 

a = 2-e- 2L2 (5.17a) 
b = c = d = e~ 2L2 (5.17b) 

In order to characterize the different phases, it is useful to compute the parameter 

A = 0-2 + f~ * = e 2 ^ 2 - 1 (5.18) 
2(ab + cd) y ' 

When < e~ 2L2 < |, we have that A > 1 and hence, the eight- vertex model is ordered 
(it belongs to the ferroelectric Regime I in Baxter's notation [3], as a > b + c + d). On 
the other hand, when e~ 2L ' 2 > |, we have that |A| < 1, and the system is disordered 
(it belongs to Baxter's Regime III). Thus, there is a critical point at the symmetric 
point fl5TT2|) . 

In order to obtain the value of the parameter /x at this critical point, we first 
need to make use of the rearrangement procedures explained in Baxter's book [21 Sec- 
tion 10.11], so that the transformed weights ( d r ) belong to Baxter's principal 
regime IV (the anti- ferroelectric phase with c > a + b + c and A < — 1). To achieve 
this, we first use the duality (or weak-graph) transformation [31 Eq. (10.2.5)]: 

a' = -(a + b + c + d) = 1 + e~ 2L2 (5.19a) 

V = ^(a + b-c-d) = l-e~ 2L2 (5.19b) 

c = ^(a-b + c-d) = l-e~ 2L2 (5.19c) 

d! = ^(a-b-c + d) = l-e" 2L2 (5.19d) 

so that the model with weights (a', b', c', d') belongs to the ferroelectric regime I (a' > 
b' + d + d' with A > 1). Then, we make the change (a r ,b r ,c r ,d r ) = (d,d',a',b'), 



39 



so that the model now belongs to the principal regime. The critical value of the 
parameter /i is given in terms of the weights (a r , b r , c r , d r ) when e~ 2Li = 1/2: 



tan^ = J**L = ^3, (5.20) 
2 V a r o r 



and hence, 

A* = f ■ (5-21) 

In this case, the critical exponents are those of the eight-vertex model for the same 
value of fi [31 Equations (10.12.24)/(10. 12.25)]: 

i/ = - (5.22a) 

y 

a = 2 - - (5.22b) 

y 

f> = 1 (5.22c) 
j = I (5.22d) 

= ^ < 5 ' 22e > 
y - ^ (5.22f) 

where y is given by (15. 5p . The numerical values for these exponents at the symmetric 
point (15.121) are: 

3 1 3 1 91 5 

v = - , a = - , (3 = - , /?' = -, 7 = — , V = - . (5.23) 

4 ' 2 ' ^ 32 ' ^ 8 ' ' 16 ' 4 v ; 

These results agree with the result found by Hintermann and Merlini [2]: the self-dual 
point (13.51) of the union-jack-lattice HMBW model is critical with critical exponent 
a = 1/2. 

We find that, on one hand, the value of fj, = 2vt/3 at the symmetric point 
e -2L 2 _ e -2i 2 = 1/2 is the same, independently of how we approach that point, 
either along the self-dual curve (15.11) when K4 — > —00, or along the HMBW-model 
line (I4.35p . However, we get two different sets of critical exponents: if we follow the 
self-dual curve (15. ip . we should use the expressions (15 .4p obtained by renormalization- 
group arguments; but along the line (I4.35p . one should use instead the eight- vertex 
exponents (15.221) . 

The exponents 1/^ = 1/2 in Eq. (15.101) and 1/v — 4/3 in Eq. (15.231) correspond to 
the renormalization exponents along the diagonal direction (i.e. along the self-dual 
curve in the large-coupling limit) in Figure [TDl and the perpendicular-to-diagonal 
direction (i.e., along the dashed line in Figure [i~2]) . respectively. The same behavior 
will be found in the triangular and hexagonal lattices (see the next two sections). 
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Finally, all points (except point B and its symmetric counterparts) on the solid 
curves in Figure [12] belong to the Ising universality class. We have checked this fact 
numerically by computing via MC simulations the values of the critical exponents on 
several points on these curves. 

5.2 Triangular- lattice ICLAT model 

The phase diagram for the symmetric AT model on the triangular lattice is qualita- 
tively different from that for the square lattice. The main reason is that the triangular 
lattice is not bipartite; therefore, there is no symmetry K2 — > —Ki- However, if we 
focus on the ferromagnetic regime K2 > 0, then it is rather similar, except that phase 
IV does not exist. See [50] for a recent study of this model. 

Another difference with respect to the same model on the square lattice is that 
the symmetric AT model on the triangular lattice is not self-dual, and it does not 
satisfy in general the star-triangle equation. However, it does satisfy the star-triangle 
equation on a certain curve in the {K21K4) plane. Temperley and Ashley [5TJ found 
that this curve is given by 

e -4^4 = 1 ( e 4K 2 _ -q (5 24) 

The qualitative behavior of the symmetric triangular-lattice AT model on this curve 
is similar to that for the square lattice: 

• The model is critical only for K4 < |log2, and it can be mapped onto the 
Gaussian model with central charge c = 1. 

• At K4 = K 2 = \ log 2, the critical curve splits into two Ising-like critical curves. 
This point corresponds to the critical coupling for a ferromagnetic 4-state Potts 
model on the triangular lattice. 

• For K4 > 1 log 2, the curve is no longer critical. 

One of the Ising-like critical curves emerging at K4 = K2 = \ log 2 goes to the 
Ising critical point at (^2,-^4) = (0, ^log3). The other line goes to the Ising critical 
point at K4 — > 00 with K 2 = \ log 3. 

As the subspace K2 = corresponds to an Ising model in the variables <tt, then 
we expect the above ferromagnetic critical point at K4 = ilog3, and an antiferro- 
magnetic critical point in the limit K4 — > —00. Thus, there is no AFIs-A curve, like 
in the phase diagram for the square lattice (see Figure fT0|) . 

Finally, the 4-state antiferromagnetic Potts model on the triangular lattice has a 
zero-temperature critical point. This point is recovered when K 2 = K4 — > — 00. 

We expect that the same critical exponents will be found for the part of the 
curve (I5.24p between the critical points for the 4-state Potts and the two decoupled 
Ising models. The effective number of Potts states along the curve (15.241) is given 
by [5U Equation (18)]: 

^ = ^ u.iV ' (5 - 25) 
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and again the relation between q and // is given by ( 15. 9p . Then, the 4-state Potts 
critical point is defined by e 4K2 = 2; hence, q = 4 and \i = 0. The decoupled Ising 
models are obtained for e 4K2 = 3, and thus, q = and \i = tt/2. These values of q, 
and /i indeed agree with those obtained for the square-lattice AT model, as we expect 
universality to hold in the ferromagnetic regime. 

However, when we take larger values of fi > tt/2, we enter the BK curve vbk for 
the triangular lattice jl9]. This is the middle branch of the curve [HI 



v 3 + 3v 2 = q, (5.26) 

and hits the v = — 1 line at g*(tri) = 2. This corresponds to y/q = —y/2, fi = 37r/4, 
and y = 3/2. These values differ from those of the square-lattice model, and lead to 
a different set of critical exponents for this point: 

v = oo, - = -2, - = -, — = -, - = t, — = 1- (5.27) 
v v 8 is 2 is 4 is 

Saleur [12] indicates that at this point fi = 37r/4 or g = 1 one has a Kosterlitz- 
Thouless transition. In Refs. [5U|55] the values for a /is, j/u, and (3/v in (I5.27P were 
obtained by studying the critical 0(2) loop model on the hexagonal lattice. It is 
interesting to note that the triangular-lattice Ising q = 2 antiferromagnet at zero 
temperature v — — 1 is critical with c = 1. 

Let us now consider the triangular-lattice ICLAT model. Its phase diagram is 
given in Figure [TBI Again it is symmetric under the interchange -H- L' 2 , so we will 
focus on the subspace L' 2 > Li- 

First of all we would like to compute the symmetric point in the (e~ 2L2 , e~ 2L2 ) 
space where the critical curve ( I5.24p ends. This is easily obtained from ( I5.24p if we 
take K±, —K 2 —> — oo with K 2 + = L 2 = L' 2 fixed. The result is 

^■ e ^) = (^^)- (5 - 28) 

This corresponds to point B in Figure [133 

Some important particular cases are given by 

• e _2L 2 = corresponds to an Ising model on the a variables. Thus, there are 
phase transitions at 

e" 2i2 = ^ (5.29a) 

e - 2L2 = +oo (5.29b) 

The first solution corresponds to the ferromagnetic Ising critical point; while 
the second one is the antiferromagnetic critical point at zero temperature. 

• The triangular-lattice HMBW model (or simply the BW model) [HOGS] cor- 
responds to the circle (I4.34p . For e~ 2Li > ^= the model is disordered, for 
e -2L 2 < moc i e i j s ordered, and finally at e~ 2L2 = the model is 
critical. 
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Figure 13: Phase diagram on the triangular-lattice ICLAT model ( 13.151) in the plane 
(e~ 2L2 , e _2L 2). Point B (depicted as a black dot) has coordinates (l/v2, l/y/2), and 
corresponds to the point where the self-dual curve of the symmetric AT model (15.24ft 
hits this plane. The ferromagnetic Ising critical point for e~ 2Li = is located at 
e -2L 2 = 1/^/3. The dashed arc going through point B corresponds to the triangular- 
lattice HMBW model, and the dashed line joining (0,1) and (1,0) corresponds to 
the bisected-hexagonal-lattice HMBW model. The notation for curves, symbols, and 
phases is as in Figure [I5J The main diagonal starting at point B (depicted as a 
blue dot-dashed line) corresponds to a line of critical points belonging to the same 
universality class as point B. 



• In Section I4.3[ we showed that the bisected-hexagonal-lattice HMBW model 
corresponds to the line (14.351) . 

Wu [20] suggested that the triangular AT model should correspond to a uniform 
8-vertex model on the kagome lattice. The details were worked out by Temperley and 
Ashley [51]. Theorem 14. II tells us that the triangular-lattice AT model with weights 
{oJk} is equivalent to a mixed AT model on the triangular and hexagonal lattices 
with weights {Qk} given by (14.81) . Then, Theorem 14.41 relates this mixed AT model 
with couplings {K2, K 2 , K^}, with an ICLAT model on the triangular lattice with 
couplings {L 2 ,L' 2 } given by (I4.23H . Notice that we have U3 = in the triangular- 
lattice AT model, so that Q2 = ^3 in the mixed AT model, which also means that the 
couplings for this model satisfy K' 2 = K4. 

We also know by following Wu [20] that this mixed AT model can be related to 
an 8-vertex model on the medial of the triangular lattice (i.e., a kagome lattice) with 
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weights 
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e K 2 +2K^ 


(5.30a) 


e -K 2 


(5.30b) 


e -K 2 


(5.30c) 




(5.30d) 



where we have taken into account that K' 2 = K4. Baxter |3] found that this uniform 
8-vertex model on the kagome lattice can be solved by transforming it to a square- 
lattice 8-vertex model when the following condition is fulfillecfl 

d c 2 -d 2 

a = (5 - 31) 

When b = c [cf. f)5.30p ]. this condition implies that ad = b 2 , or equivalently, K2 = 0. 
Therefore, the 8-vertex model on the kagome lattice is soluble by Baxter's methods 
if the weights have the form: 

a = e 2j? 2 (5.32a) 
b = c = 1 (5.32b) 
d = e~ 2R i (5.32c) 
This 8-vertex model is critical if a = b + c + d. This equation implies that 

sinh(2^) = 1 & e- 2L2 = e' 2 ^ = — . (5.33) 

This corresponds to the symmetric point (15.281) . 

To obtain the value of \i corresponding to this point, we proceed as for the square 
lattice. If e~ 2L2 > l/y/2, the model is disordered (it belongs to regime III). Using the 
weak-graph transformation ( 15. 19ft . we can get a point in regime I with weights: 

a' = 1 + cosh 2K' 2 (5.34a) 
b' = c ' = sinh2^2 (5.34b) 
d! = -1 + cosh 2K' 2 (5.34c) 



6 Notice we have followed Wu's notation 20 to go from spin configurations to arrow configu- 
rations, and this choice is non-equivalent to that of Baxter [3] or Temperley-Ashley [ST]. We can 
"translate" Baxter's formulas to our notation by using the transformation (a, &, c, d) — > (c, d, a, b). 
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If we do the transformation (a', b', c', d!) — > (c', d 1 , a 1 , b') = (a r , b r , c r , d r ), we get finally 
into the principal regime, with weights 

a r = d r = sinh2i^2 

b r = -1 + cosh 2K' 2 

c r = 1 + cosh 2K' 2 

Using these weights, we compute the critical value of fi: 

tanf = = V2 + 1, 

2 V a rb r 

and hence, 

3tt 

" = T" 

The numerical values for these exponents at the symmetric point ( I5.28P are: 

" = i- ^Tr r-h- -< = l- y = 5- (5 ' 38) 

These results agree with the result found by Baxter and Wu [4"ll5ll56]. 

It is worth noticing that these exponents coincide with those for the ferromag- 
netic 4-state Potts model \i = 0. In fact, these two models are believed to belong 
to the same universality class [6]. One striking difference between the BW model 
and the 4-state Potts model is that the latter displays logarithmic corrections (both 
multiplicative and additive) [814TO]. while the former does not have any of these. 

A qualitative phase diagram for this model is shown in Figure [T3J We expect, 
as in the square-lattice case, that at the point B (where the self-dual curve of the 
symmetric AT model hits the ICLAT plane), the critical curve splits into two curves 
going to Ising-like critical points on the axis. These curves are expected to be Ising- 
like. We also know that the line (14.351) corresponds to the HMBW model on the 
bisected-hexagonal lattice. The phase diagram shows that this model has two critical 
points, contrary to what we know about the same model on the other two triangula- 
tions (triangular and union-jack lattices). Our numerical MC results show that their 
location iso 

e 2L 2 = o.4196(15), e 2L 2 = 0.5803(8) (5.39a) 
e 2L2 = 0.5803(8) , e 2L ' 2 = 0.4196(15) (5.39b) 

In terms of the HMBW coupling J (I4.33p . we obtain the critical points 

J lc = 0.5003(5), J 2 , c = 0.3857(4). (5.40) 

Indeed, we can check that these two values are mutually dual: using the definition 
f)3.4p . we arrive at vxv 2 = 2.001(3). Therefore, we suggest the following: 

7 We performed the MC computation on the kagome lattice (see section HT4|) . The values of the 
couplings for the triangular lattice are the square of those obtained for the kagome lattice. 



(5.35a) 
(5.35b) 
(5.35c) 

(5.36) 
(5.37) 
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Conjecture 5.1 The Hintermann-Merlini-Baxter-Wu model A3. l)) // TOj) on the bisected- 
hexagonal lattice has two critical points 

Ji,c =\, J ^=\ lo S \JZj) W °- 3859684164 • ( 5 - 41 ) 

These two points are mutually dual. 

Finally, the diagonal line starting at point B and depicted as a dot-dashed line 
on Figure [[3] belongs to the same universality class as point B. This is supported 
by our MC simulations along this line. The above statement is true for all points 
on this diagonal line, except in the limit e L2 = e L ' 2 — > +00, where we recover a 
hexagonal-lattice FPL model with fugacity n = 2. This FPL model is equivalent 
to the zero-temperature limit of the triangular-lattice 4-state Potts antiferromagnet 
[STJEB]. Therefore, its continuum limit can be described by a conformal field theory 
of central charge c = 2 [25J. 

5.3 Hexagonal-lattice ICLAT model 

The phase diagram for the symmetric AT model on the hexagonal lattice is quali- 
tatively the same as for the square lattice, as both lattices are bipartite. However, the 
AT model on the hexagonal lattice is not self-dual; but it does satisfy the star-triangle 
equation on a certain curve in the (K 2 , K4) plane. Temperley and Ashley [51] found 
that this curve is given by 

e~ 4Ki + e~ 2Ki cosh 2K 2 = sinh 2 2K 2 (5.42) 

We then find that 

• The model is critical only for K\ < -log 5, and it can be mapped onto the 
Gaussian model with central charge c = 1. 

• At K4 = K 2 = ^log5, the critical curve splits into two Ising-like critical curves. 
This point corresponds to the critical coupling for a ferromagnetic 4-state Potts 
model on the hexagonal lattice. 

• For K4 > I log 5, the curve is no longer critical. 

As the subspace K 2 = corresponds to an Ising model in the variables err, then 
we expect the above ferromagnetic critical point at K A = |log(2 + a/3), and an 
antiferromagnetic critical point at — — | log (2 + v3)- There is a curve of Ising-like 
critical points emerging from this point, and going to K4 — > —00, as for the square 
lattice. 

Finally, the 4-state antiferromagnetic Potts model on the hexagonal lattice is 
disordered at all temperatures [SSJED], thus none of the above critical curves cross 
the line K 2 = -K 4 > 0. 

We expect that the same critical exponents will be found for the part of the curve 
(15.421) between the critical points for the 4-state Potts and the two decoupled Ising 
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models. The effective number of Potts states along the curve (I5.42p could be obtained 
from the results given in Ref. [51J: 

y/q = y/l + OKij [£(K 2 ) - 2] , (5.43) 

where £(AV) is given by 



4tanh i (2A 2 ) 



(5.44) 



2tanh 2 (2A 2 ) 

Again the relation between q and \x is given by ( 15. 9p . Then, the 4-state Potts critical 
point is defined by e 2K2 = a/5; hence, £ = 3, q = 4 and fj, — 0. The decoupled Ising 
models are obtained for e 2K2 = 2 + \/3, and thus, £ = 2, g = and [i = ir/2. These 
values of q, and /i indeed agree with those obtained for the square-lattice AT model, 
as we expect universality to hold in the ferromagnetic regime. 

When we take larger values of fi > n/2, we enter the BK curve v-bk for the 
hexagonal lattice. This is the lower branch of the curve [4"5 | l52 | l53] 

v 3 -3vq = q 2 , (5.45) 

and hits the v — — 1 line at q*(hex) = (3 — \fh)/2. This corresponds to y/q = 
(1 — \/5)/2, /i = 37r/5, and y = 6/5. This value coincides with the one we obtain 
when we take the limit K 2 — > +oo in (I5.43p . as lim^^oo £(A" 2 ) = (1 + \/5)/2. These 
values lead to the following set of critical exponents for this point: 

V = V a = ~l P = I ' = 7 = V 7 ' = T' (5 ' 46) 

Let us now consider the hexagonal-lattice ICLAT model. Its phase diagram is 
given in Figure [HJ Again it is symmetric under the interchange L 2 -H- L' 2 , so we will 
focus on the subspace L' 2 > L 2 . 

First of all we would like to compute the symmetric point in the (e~ 2L2 , e~ 2L ' 2 ) 
space where the critical curve ( I5.42p ends. This is easily obtained from ( I5.42p if we 
take A4, — K 2 — > —00 with A 2 + K4 = L 2 = L' 2 fixed. The result is 

(e-,e-) = (^.^) ■ (5.47) 

This corresponds to point B in Figure HU and it is equivalent to point B' of coordinates 

e- 2L2 ,e~ 2L '^ = (I + a/5,1), (5.48) 

because of the symmetry (I3.22p . 

Some important particular cases are given by 

• e~ 2L ' 2 = corresponds to an Ising model on the a variables. Thus, there are 
phase transitions at 

e -2L 2 = 2 ±y/3. (5.49) 

The '+' solution corresponds to the ferromagnetic Ising critical point; while the 
'— ' solution is the antiferromagnetic critical point. Symmetry f !3.22|) implies 
that these points are mapped to points at infinity with slopes 2 ± y/3. 
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Hex-lattice ICLAT model 




Figure 14: Qualitative phase diagram on the hexagonal-lattice ICLAT model (13.151) in 
the plane (e~ 2L2 , e~ 2Ij2 ). Point B of coordinates ((v5 — 1)/4, (y/E— l)/4) corresponds 
to the point where the self-dual curve of the symmetric AT model C I5 .421) hits this 
plane, and the other two points are equivalent to the former one because of the 
symmetries of this model. The ferromagnetic and antiferromagnetic Ising critical 
points for e~ 2Ij2 = are located at 2 ± a/3, respectively. The dashed arc going 
through point B corresponds to the models that are solvable using Baxter's method. 
The dashed line corresponds to the uniform HMBW model on the bisected-hexagonal 
lattice. The notation for curves, symbols, and phases is as in Figure [12j 



• The point (e~ 2i2 , e~ 2L ' 2 ) = (1,1) corresponds to the zero-temperature 4-state 
antiferromagnetic Potts model on the hexagonal lattice. This model is known 
to be disordered [59|l60]. 

• In Section 14.31 we showed that the bisected-hexagonal-lattice HMBW model 
corresponds to the line (I4.35p . In agreement with Conjecture 15.11 we find two 
critical points on this line at the same values as for the triangular-lattice ICLAT 
model (I5.39P (recall that in Section 14.31 we obtain maps from the bisected- 
hexagonal-lattice HMBW model to the ICLAT model on decorated triangular 
or hexagonal lattices). 

The kagome lattice is also the medial of the hexagonal, thus we can repeat the same 
argument to relate the hexagonal-lattice ICLAT model to the kagome-lattice 8-vertex 
model [5T], using a mixed AT model as the intermediate step (20] • Again, Theorem l4.ll 
tells us that the hexagonal-lattice AT model with weights {u^} is equivalent to a 
mixed AT model on the hexagonal and triangular lattices with weights {Qk} [cf. 
(I4.8p ], and this latter model us related to a hexagonal-lattice ICLAT model with 
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couplings {L2,L' 2 } [cf. f!4.23p ]. Notice that we have U3 = in the hexagonal-lattice 
AT model, so that Q 2 — ^3 i n the mixed AT model, which also means that the 
couplings for this model satisfy K 2 = K A . 

The results of Wu [20J ensure that this mixed AT model can be related to an 
8-vertex model on the medial of the hexagonal lattice (i.e., a kagome lattice) with 
weights given by (I5.30p . 

The only difference with respect to the triangular-lattice case is that in this latter 
case, the a spins live on the triangular lattice, and the r spins on the dual hexagonal 
lattice; while for the hexagonal-lattice case, the situation is reversed: the a spins live 
on the hexagonal lattice, and the r spins live on the dual triangular lattice. This 
implies that the condition for the corresponding kagome-lattice 8-vertex model to be 
soluble by Baxter's method now readstJ 

c d 2 - c 2 

(5.50) 



a a 2 — b 2 

When b = c [cf. (15.30!) ]. this condition implies that the couplings K 2 and K' 2 should 
satisfy the following equation: 

x 3 e 2K> _ x 2 e AK> _ x JK' _|_ ]_ = Q , (5.51) 

where 

x = e~ 2K2 . (5.52) 

It is difficult to write the solutions of the above equation (I5.5ip as simple expressions 
in the ICLAT couplings L 2 and L' 2 . However, it is not difficult to prove that the 
following points provide solutions of (|5.51jl : 



-*,e->) = P^,^J (5.53b) 

Therefore, the symmetric point (15. 47ft belongs to this curve where the ICLAT model 
is soluble 8-vertex model on the square lattice. 

Furthermore, the symmetric point (I5.47P is the only one where the 8-vertex model 
is critical. The condition for this 8-vertex model to be critical is 

a = b + c + d, (5.54) 
which in terms of the couplings K 2 , K' 2 means that 

sinh2^2 = e" 2 ^ 2 . (5.55) 



8 Again, we have followed Wu's notation [20] to go from spin configurations to arrow con- 
figurations. We can "translate" Baxter's formulas to our notation by using the transformation 
(a, b, c, d) — > (d, c, a, b). 
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This condition implies that 




e - 2L 2 = tanh2K 2 = e~ 2L2 . (5.56) 

Therefore, the model is critical precisely at the symmetric point (I5.47p . 

To obtain the value of \x corresponding to this point, we proceed as for the trian- 
gular lattice. By rearranging the weights appropriately so that the model is mapped 
to the principal regime, we obtain for the value of fi at the critical point (I5.47P 

tan | = \l - ' ; v " , (5.57) 
and hence, 

A* = y- (5-58) 
The numerical values for these exponents at the symmetric point ( I5.47P are: 

"=-, « = -, P = —, P = ~, 7=-, i = -■ (5.59) 
6 3 48 6 ' ' 24 ' ' 3 V ; 

Again we find that this set of critical exponents do not agree with the set (I5.46P 
obtained by following the critical curve of the self-dual curve of the hexagonal-lattice 
symmetric AT model. 

Finally, as for the square lattice, we find numerically in our MC simulations that all 
solid curves in Figure O are Ising-like, except point B and its symmetric counterparts. 



5.4 Kagome- lattice ICLAT model 

The phase diagram for the symmetric AT model on the kagome lattice is expected 
to be qualitatively similar to that of the same model on the triangular lattice. The 
main reason is that both lattices are not bipartite. Unfortunately, the exact location 
of any of the critical curves is unknown for this case. See Ref. (50] for a recent study. 

What is exactly known about the behavior of the symmetric kagome-lattice AT 
model on this curve reduces to the solution of the Ising model |61J : 

• For K 2 = 0, it is critical in the ferromagnetic regime at coupling 

K 4 = -log(3 + 2V3). (5.60) 

• In the antiferromagnetic regime it is always disordered. 

Finally, the 4-state Potts model on the kagome lattice has a ferromagnetic crit- 
ical point at K 2 = K± ~ 0.2873, and in the antiferromagnetic regime, it is always 
disordered. 

Let us now consider the kagome-lattice ICLAT model. Its phase diagram is given 
in Figure [T5j Again it is symmetric under the interchange L 2 L' 2 , so we will focus 
on the subspace L' 2 > L 2 . 
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The position of the critical point on the symmetric subspace is not exactly known 
for this model. Our numerical MC results show that this point has coordinates 
(0.45411(3), 0.45411(3)). This corresponds to point B in Figure [TBI 

Some important particular cases are given by 

• e~ 2L 2 = corresponds to an Ising model on the a variables. Thus, there are 
phase transitions at 

e -2L 2 = 1 _ , 561) 

V3+2V3 

• In Section I4.3[ we showed that the uniform bisected-hexagonal-lattice HMBW 
model corresponds to the circle f)4.34p . In agreement with Conjecture 15.1} we 
find that the critical curve for the kagome-lattice ICLAT model crosses the 
above circle at two points: 

e 2La = 0.6478(5), e 2L * = 0.7618(3) (5.62a) 
e 2i2 = 0.7618(3), e 2L * = 0.6478(5) (5.62b) 

Indeed, the values of the HMBW coupling J are the same as in (I5.40p . 




Figure 15: Phase diagram on the kagome-lattice ICLAT model (13.15!) in the plane 
(e~ 2L2 ,e~ 2L 2). Point B has coordinates (0.45411(3), 0.45411(3)), and corresponds to 
the point where the self-dual curve of the symmetric AT model hits this plane. The 
ferromagnetic Ising critical point for e~ 2Ij2 = is located at e~ 2L2 = 0.39360791. 
The dashed arc corresponds to the bisected-hexagonal-lattice HMBW model. The 
notation for curves, symbols, and phases is as in Figure [£21 
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6 Discussion 



The main goal of this paper is the study of the infinite-coupling-limit Ashkin- 
Teller model on several common lattices, namely, square, triangular, hexagonal, and 
kagome. We also have considered a generalization of the well-known Baxter-Wu 
model: the Hintermann-Merlini-Baxter-Wu model defined on any plane Eulerian 
triangulation. 

We have first investigated the relations between these two models. As a side effect 
we have also considered the Ashkin-Teller and the mixed Ashkin-Teller models, which 
are known in the literature. In particular we find that 

• The partition functions for the AT model on a graph G ( 13. 6 p and the mixed 
AT model (I3.25P on G and its dual G* are equal (modulo some unimportant 
multiplicative factors): 

Zat(G] LUq, Ui, L} 2 , U} 3 ) (X Z m AT(G ! , G*] CUo, U>1, U)2, u 3 ) , 

where the weights Qk are given in terms of the weights by ( 14. 8p . (See The- 
orem STTJ) As a side effect, we recover the duality relation for the AT model 
(jUgp . (See Theorem Ed) 

• The partition function for the ICLAT model on a graph G ( I3.16P and the par- 
tition function of the mixed AT model (I3.25P on G and G* are again equal 
(modulo some unimportant multiplicative factors): 

Z m AT{G, G*; Cc>o, Wi, ui 2} uj 2 ) oc Ziclat(G; L 2 , L' 2 ) 

where the weights couplings L 2 ,L 2 are given in terms of the couplings K 2 , K 2 
(or the weights Qk) by ( I4.23j) . (See Theorem 14.31 ) 

• Using the partial-trace transformation (Lemma 14. 5[) we obtain that the HMBW 
model on certain lattices is a subspace of the ICLAT on other lattices. In 
particular Corollary 14.71 implies: 

— The triangular-lattice HMBW is equivalent to a triangular-lattice ICLAT 
model on the subspace (14 .34 p . 

— The bisected-hexagonal-lattice HMBW is equivalent to a kagome-lattice 
ICLAT model on the subspace (I4.34|) . 

— The union-jack-lattice HMBW is equivalent to a square-lattice ICLAT 
model on the subspace (14 .35 p . 

— The bisected-hexagonal-lattice HMBW is equivalent to a hexagonal-lattice 
ICLAT model on the subspace (I4.35p . 

— The bisected-hexagonal-lattice HMBW is equivalent to a triangular-lattice 
ICLAT model on the subspace (I4.35p . 
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• We have found in Section FOl that the partition function for the BW model (i.e., 
the triangular-lattice HMBW model (I3.2p ) can be written as a partition function 
of a 2-color-non-overlapping Eulerian bond model on a hexagonal lattice with 
weights 

The second important result is the computation of the phase diagram of the 
ICLAT model on the square, triangular, hexagonal, and kagome lattices. In the first 
three cases, we know the exact point (labeled B in the figures) where the self-dual 
curve of the corresponding symmetric AT model hits the ICLAT model in the limit 
of infinite couplings. Using the expression for the critical exponents of the symmetric 
AT along this self-dual curve, we are able to make predictions about the critical 
exponents at B. On the other hand, in the ICLAT plane, we can also use the relation 
between the AT model (and its infinite-coupling limit) and the 8-vertex model. From 
this relation another set of critical exponents arise. In the case of the square and 
triangular lattice, these predictions agree with the results of Hintermann-Merlini [2] 
and Baxter- Wu [HE]. It is interesting to note that these two sets of critical exponents 
do not coincide. 

Finally, the exact location of the self-dual curve of the kagome-lattice AT model 
is not known. In this case, we have obtained a high-precision Monte Carlo estimate 
for the positions of point B and of the two points where the critical curves cross 
the bisected-hexagonal HMBW subspace. From these latter values, we have deter- 
mined the two critical points of the bisected-hexagonal-lattice HMBW model (see 
Conjecture 15. ip . 

A Decimation transformation 



01 0"2 01 cr 2 




Figure 16: (a) The HMBW model with generic 3-spin interaction Ji0o0i02 + 
J20~oo~20~3 + </30oO"3CX4 + J 4000 40i ■ (b) The AT model with 2-spin interaction if 2 0204 + 
K^aias, and 4-spin interaction K^ai ^u^a^ obtained from (a) after summing over the 
values of 00 = ±1. 

In this appendix we will consider the decimation transformation, that can be 
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considered as a generalization of the k = 2 case of the partial-trace transformation 
dealt with in Section 14.31 

Let us suppose we have a general (i.e., face- dependent coupling) HMBW model 
defined on an Eulerian plane triangulation G = (V,E). Then the vertex set has a 
tripartition V — V\ U V 2 U V3. Let us further suppose that the vertices belonging to 
the sublattice V 3 have all degree A3 = 4. Then, this general HMBM model can be 
transformed into a mixed AT model by performing the sum over all spins living on 
V 3 . (See Figure [TBI) 

2J exp [cr ( Ji(Xi<72 + J 2 o 2 o 3 + + Ja<?±<?i)\ 

<x =±l 

= 2 exp [K + K 2 a 2 a 4: + K' 2 Oi(j 3 + K A aia 2 a 3 a^ . (A.l) 

Inspecting all the possible configurations for the spins cr^ with i = 1,2,3,4, one has 
the following 4 independent equations 

coo = exp(K + K 2 + K' 2 + K A ) = cosh(Ji + J 2 + J 3 + J 4 ) (A.2a) 

wi = exp(K - K 2 + K' 2 - K A ) = cosh( J x + J 2 - J 3 - J 4 ) (A. 2b) 

u 2 = exp(K + K 2 -K' 2 - K 4 ) = cosh( J x - J 2 - J 3 + J 4 ) (A.2c) 

cu 3 = exp(A - K 2 - K' 2 + A 4 ) = cosh( Ji - J 2 + J 3 - J 4 ) (A. 2d) 

These equations have the solution: 

e 4Ko = cosh( Ji + J 2 + J 3 + J 4 ) cosh( J x + J 2 - J 3 - J 4 ) 

x cosh( — J 2 — -^3 + Ja) cosh( — J 2 + J 3 — J 4 ) (A. 3a) 

e 4K 2 _ cosh( Jj + J 2 + J 3 + J 4 ) cosh( Ji -J2-J3 + Ja) (A 3b) 

cosh( Ji + J 2 — J 3 — J 4 ) cosh( Ji — J 2 + J 3 — J 4 ) 

g4K , _ cosh( J x + J 2 + J 3 + J 4 ) cosh( Ji + J 2 - J 3 - J 4 ) (A 3c) 

cosh( Ji — J 2 — ^3 + Ja) cosh( Ji — J 2 + J 3 — J 4 ) 

g4 ii- 4 _ cosh( Jx + J 2 + J 3 + J 4 ) cosh( Ji - J 2 + J 3 - J 4 ) (A 3d) 

cosh( Ji + J 2 — J 3 — J 4 ) cosh( J x — J 2 — J 3 + J 4 ) 

Notice that the 2-spin terms couple spins placed on the diagonal; not nearest-neighbor 
ones. Therefore, we obtain a mixed AT model on the graph G\ = (Vi,Ei) [defined 
by the vertices in Vi], and its dual G\ = (V 2 , El) [defined by the vertices in V 2 \- 

Z nMBW (G;J) = 2\ V ^Z^ T {G X ,G\-K 2 ,K' 2 ,K A ) (A.4a) 
= ^at^G*; u^,^, u, 3 ) (A.4b) 

where the relation between the mixed AT couplings and the original HMBW cou- 
plings is given in (1A.3|) . and the weights u}^ are given in (IA.2I) . To emphasize that 



the couplings are face- or edge-dependent, we use boldface letters. The relation be- 
tween general HMBW, AT, and mixed AT models can be generalized easily from the 
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mappings discussed in Section HI Indeed, we can interchange V\ and V 2 in the above 
discussion. 

We can rewrite (lA.4j) in terms of the partition function of an AT model on G\ by 
using (flTTUjl : 

Z UMBW (G;J) = 2\ v ^\ +l Z AT {G\- Q ,Q 2 ,Q U u> 3 ) , (A.5) 

where the new weights are given by 

cD = ujq + u 2 = 2 cosh( J 4 + J 2 ) cosh( J 3 + J 4 ) (A. 6a) 

cDi = ujq — u 2 = 2 sinh( J\ + J 2 ) sinh( J 3 + J 4 ) (A. 6b) 

0)2 = wi + W3 = 2cosh(Ji — J 2 ) cosh(J 3 — J 4 ) (A. 6c) 

0)3 = Wi + = 2sinh(Ji — J 2 ) sinh(J 3 — J 4 ) (A.6d) 

The most interesting case corresponds to a homogeneous HMBW model (13.1 p / (13.21) 
with coupling J. Then, (1A.5|) reduces to a ICLAT model: 

Z UMBW (G;J) = 2l y3 l-l y2 l +1 Z AT (Gt;2cosh 2 (2J),2sinh 2 (2J),2,0) (A.7a) 
= Z lCLAT (Gl;L 2 ,L' 2 ) (A.7b) 

where the couplings L 2 , L' 2 are given by: 

e~ 2L2 = (A.8a) 

cosh (2 J) V ; 

e" 2L 2 = tanh 2 (2J) (A.8b) 

Indeed, the uniform HMBW model on an Eulerian planar triangulation G maps into 
the line ffl~35l) of the ICLAT model on the graph G* x . 

If we apply this transformation to the union-jack lattice we get case (b) of Sec- 
tion I4.3[ while if we apply it to the bisected-hexagonal lattice, we obtain cases (c) 
and (d) of Section 14.31 
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